AD“A207  549 


inc  ffiE  mp» 


ESMC-TR-8a-03 


A  REVIEW  OF 

ENERGY  RELEASE  PROCESSES 
FROM  THE  FAILURE  OF 
PNEUMAnC  PRESSURE  VESSELS 


August  1988 


s 


DTIC 

ELECTS 
MAY  0  1 1989 


(ijo  H 


Geoeral  Physics  CorporatioD 
Cape  CAnareral  Air  Force  Station 
Florida 


Eastern  Space  and  Missile  Center 
Air  Fbrce  Systems  Command 
Fitrick  Air  Force  Base,  Florida 


APPROVED  FOR  PUBLIC  RELEASE 
DISTRIBUTION  UNLIMITED 


89  5  01  o85 


TfflS  DOCUMENT  IS  BEST 
QUALITY  AVAILABLE.  THE  COPY 
FURNISHED  TO  DTIC  CONTAINED 
A  SIGNIFICANT  NUMBER  OF 
PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


NOTICE 


When  governaent  specifications,  or  other  data  are  used  for  any  purpose  other 
than  in  connection  with  a  definitely  related  gOTernaent  procureaent  operation, 
the  United  States  Govemaent  thereby  incurs  no  responsibility  nor  any  obligation 
whatsoever;  and  rha  fact  that  the  goTemaent  aay  have  fomilated,  furnished,  or 
in  any  way  supplied  the  said  specifications,  or  other  data,  la  not  to  be 
regarded  by  iaplication  or  otherwise  as  in  any  aanner  licensing  the  holder  or 
any  other  person  or  corporation,  or  conveying  any  right  or  permission  to 
manufacture  use,  or  sell  any  patented  invention  that  aay  in  any  way  be  related 
thereto. 

This  report  has  been  reviewed  by  the  Office  of  Public  Affairs  (ESMC/PA)  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS).  At  NTIS,  it 
will  be  available  to  the  general  public.  Including  foreign  nations. 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


£i^ 


Project  Engineer 


FOR  THE  DIRECTOR: 


COL,  USAF 

Director  of  Safety 


THOMAS  L.  WILSON,  LTC,  USAF 
Chief,  Missile  Systems  Safety 
Division 


«ieuW>TV~e~L*$«l»'C»tlON  0>  THI«  >»0t 


y  WI^ORT  S«CU«ITV  CUA»*<HeATION 

iiMri  /kccTPTrn 


I*.  SeCUAlTV  CLASSinCATION  authority 


3b.  oECLASSI^ICATION/OOWNanAOINaSCHtOULE 


REPORT  DOCUMENTATION  PAGE  _ _ 


lb.  HESTHICTIve  MARKINGS 

NONE  _  - 


3.  OISTHIBUTION/AVAILABILITV  of  REPORT 

APPROVED  FOR  PUBLIC  RELEASE: 
DISTRIBUTION  UNLIMITED 


4  performing  ORGANIZATION  REPORT  NUMBERISI 

GP-R-213083  _ 


Sa  name  of  PERFORMING  ORGANIZATION  Sb.  OFFICE  SYMBOL 

GENERAL  PHYSICS  CORP.  ufoppuccuic 

coarr  rnfl'tT  nPPirF 


6c.  address  Idly.  Slal*  tnd  ZIP  Codtl 

P.O.  BOX  21265  ^ 

KENNEDY  SPACE  CENTER,  FL  32815 


a*.  NAME  OF  FUNOING/SPONSORING 

organization 

ESMC 


ec  AOORESS  iCily.  Stale  and  ZIP  Codel 


Bb.  OFFICE  SYMBOL 
Ilf  applleablel 

SEM 


5.  monitoring  ORGANIZATION  REPORT  NUMBERISI 

ESMC-TR-88-03 


7».  NAME  OF  monitoring  ORGANIZATION 

RECTORATE  OF  SAFETY _ 


7b.  AOORESS  rcily.  Stale  and  ZIP  Codel 

EASTERN  SPACE  AND  MISSILE  CENTER 
PATRICK  AIR  FORCE  BASE,  FL  32925 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

F08606-86-C-0030  _ 


10.  SOURCE  OF  FUNDING  NOS. 


program 
element  no. 


PROJECT 

TASK 

NO. 

NO. 

WORK  UNIT 
NO. 


II.  title  ilncluJ*  Sccurily  CNuiA'rallonJ 

Review  of  Energy  Release  Processes  from  the 


13.  personal  AUTHORISI  ^  _ _ 

.COLEMAN,  M.CAIN,  R.OANNA,  C. HARLEY  AND  D. SHARP,  GENERAL  PHYSICS  CORPORATION 


134.  TYPE  OF  REPORT 

INITIAL  RELEASE 


16.  supplementary  notation 


13b.  time  COVERED 

FROMpGinm  toRRQ831 


14.  DATE  OF  REPORT  lYr,  Mo..  Dayl 

880831 


is.  page  count 

125 


COSATI  COOES 


18,  SUBJECT  TERMS  iContinut  on  rtxftrge  if  ncctstary  end  identify  by  num^rj 

PRESSURE  VESSELS,  BLAST  WAVE,  FR-AGMENTATION,  VESSEL 
BURST,  BURST  TEST,  AND  FRAGMENTS 


19.  abstract  iContinut  on  revoree  if  neeeuory  and  identify  by  block  numberf 

This  report  presents  the  results  of  a  review  of  the  curr^ently  available  methodologies  for 
predicting  the  characteristics  of  the  blast  wave  and  fragments  produced  by  the  failure  of 
a  pneumatic  pressure  vessel.  Methodologies  are  presented  and  evaluated  for  their  usefulness 
in  accurately  predicting  the  hazards  presented  by  a  pressure  vessel  failure.  Deficiencies 
in  the  current  methodologies  are  identified  as  well  as  the  areas  which  require  additional 
testing  for  verification.  Computer  codes  which  predict  blast  wave  overpressure  and  fragment 
initial  velocity  and  range  are  described  and  their  results  are  compared  to  experimental 
data.  Finally,  a  preliminary  test  program  is  presented  to  provide  guidance  on  setting  up 
burst  tests  to  obtain  necessary  additional  data. 


70  OiSTRIBUTION/AVAILABILITV  OF  abstract 
UNCLASSIFiEO/UNLIMITEO  (35  SAME  AS  RPT.  □  OTIC  USERS  □ 


77,  NAME  OF  RESPONSIBLE  INDIVIDUAL 


Bobby  L.  Webb 


DO  FORM  1473,  83  APR 


21.  ABSTRACT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


22b.  TELEPHONE  NUMBER  22e.  OFFICE  SYMBOL 

tfnjetude  Area  Codel 

407-494-7077  SEM 


EDITION  OF  1  JAN  73  IS  OBSOLETE,  UNCLASSIFIED _ 

security  CLASSIFICATION  OF  THIS  PAGE 


UNCLASSIFIED 


SecuBiTV  CLASSIFICATION  Of  THIS  FACE 


11.  TITLE  (Continued) 


Failure  of  Pneumatic  Pressure  Vessels  (UNCLASSIFIED) 


V  > 


Aooesslon  For 


NTIS  0RA4I 
DTIC  TAB  n 

Jui;',  ir  I  '  l- !  on - 


L-l  ‘  in/ 

A  '  ;  .>  ;  i  1  ly 

aiL;t.  fv  -ol/jl 


UNCLASSIFIED 


SECOAlTV  classification  OF  TmiS  »A0£ 


MCOKITV  CCASSH'ICATION  OF  THIS  PAGE _ 


REPC 


R6P0«T  SECUKITV  CLASSIFICATION 

UNCLASSIFIED 


3*.  SeCL'»  CLASSIFICATION  AUTHOflITY 


ab.  OECLASSIFICATION/OOWNGRAOINO  SCHEDULE 


4.  PERFORMING  ORGAFQ^i^lON  REPORT  NUMBER(S) 


REPORT  DOCUMENTATION  PAGE 


PERFORMING  ORGAFQ^i^ll 

GP-R-213083  y  V 
\r  1^ 


Ry 


6a.  NAME  OF  PERFORMING  ORGANIZATION 

GENERAL  PHYSICS  CORP. 

SPACE  COAST  OFFICE _ 


6c.  ADDRESS  (Cily.  Slate  and  ZIP  Code} 

P.O.  BOX  21265 

KENNEDY  SPACE  CENTER,  FL  32815 


b.  OFFICE  SYMBOL 
tit  applicable) 


8a.  NAME  OF  FUNOING/SPONSORING 
ORGANIZATION 

ESMC 


Sc.  ADDRESS  (Cily.  Slate  and  ZIP  Code) 


8b.  OFFICE  SYMBOL 
(If  applicable) 


1b.  RESTRICTIVE  MARKINGS 

NONE 


3.  DISTRIBUTION/AVAILABILITY  OF  REPORT 

APPROVED  FOR  PUBLIC  RELEASE: 
DISTRIBUTION  UNLIMITED 


6.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

ESMC-TR-88-03 


7a.  NAME  OF  MONITORING  ORGANIZATION 

DIRECTORATE  OF  SAFETY 


7b.  ADDRESS  (City,  Stale  and  ZIP  Code) 

EASTERN  SPACE  AND  MISSILE  CENTER 
PATRICK  AIR  FORCE  BASE,  FL  32925 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

F08606-86-C-0030 


10.  SOURCE  OF  FUNDING  NOS. 


EASTERN  SPACE  AND  MISSILE  CENTER 

PATRICK  AIR  FORCE  BASE,  FL  32925 

PROGRAM 
ELEMENT  NO. 

78022F 

PROJECT 

NO. 

TASK 

NO. 

WORK  UNIT 
NO. 

11.  TITLE  (Include  Security  Ctauificalion) 

5 

Review  of  Energy  Release  Processes  from  the 

13.  PERSONAL  AUTHOR(S) 

M. COLEMAN,  M.CAIN,  R.  DANNA,  C. HARLEY  AND  D. 


13a.  TYPE  OF  REPORT 

INITIAL  RELEASE 


13b.  TIME  COVERED 

FROM  861001 


SHARP,  GENERAL  PHYSICS  CORPORATION 


14.  DATE  OF  REPORT  (Yr.,  Mo..  Day)  15.  PAGE  COUNT 

880831  125 


17. 


FIELD 


COSATI  COOES 


GROUP 


18.  SURjECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  Woclr  number) 

^PRESSURE  VESSELS,'  BLAST  WAVE,'  FRAGMENTATION,'  VESSEL 
BURST,'  BURST  TEST,'-A5D  FRAGMENTS.  I'l.J).' - 


19.  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

^This  report  presents  the  results  of  a  review  of  the  currently  available  methodologies  for 
predicting  the  characteristics  of  the  blast  wave  and  fragments  produced  by  the  failLure  of 
a  pneumatic  pressure  vessel.  Methodologies  are  presented  and  evaluated  for  their 
usefulness  in  accurately  predicting  the  hazards  presented  by  a  pressure  vessel  failure. 
Deficiencies  in  the  current  methodologies  are  identified  as  well  as  the  areas  which 
require  additional  testing  for  verification.  Computer  codes,  which  predict  blast  wave 
overpressure  and  fragment  initial  velocity  and  range  are  described  and  their  results  are 
compared  to  experimental  data.  Finally,  a  preliminary  test  program  is  presented  to 
provide  pldance  on  setting  up  burst  tests  to  obtain  necessary  additional  data. 


20.  OISTRIBUTION/AVAILABILITY  OF  ABSTRACT  21.  ABSTRACT  SECL 

UNCLASSIFIED/UNLIMITEO  n  SAME  AS  RPT.  □  OTIC  USERS  □  UNCLASSIFIED 


22,.  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Bobby  h.  Webb 


DD  FORM  1473,  83  APR  EDITION  OF  1  JAN  73  IS  OBSOLETE. 


21.  ABSTRACT  SECURITY  CLASSIFICATION 


22b.  TELEPHONE  NUMBER 
ffnciude  Area  Code) 

407-494-7077 


22c.  OFFICE  SYMBOL 

SEM 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


UNCLASSIFIED 


96CU«lTV  CLASSIFICATION  OF  THIS  PAOf _ 

11.  TITLE  (Continued) 

Failure  of  Pneumatic  Pressure  vessels  (UNCLASSIFIED) 


UNCLASSIFIED _ 

SeCOfliTV  CLASSIFICATION  OF  THIS  FACE 


COMKEMT  FORM 


TXTUB:  R«vi«w  of  Energy  Release  Processes  from  the 
Failure  of  Pnetunatic  Pressure  Vessels 


Fublioation:  B8MC-TR-B6-03  Revision:  0  (Initial  Release) 


1.  OSAF  solicits  your  comments  concerning  this  report  so  that  its 
usefulness  may  be  liQ>roved  In  later  editions.  Send  any  comments  to 
the  following  address: 


Directorate  of  Safety 
Attn:  B.L.  Webb 

Eastern  Space  and  Missile  Center  (AFSC) 
Patrick  Air  Force  Base,  Florida  32925 


2.  Conuaents  are  solicited  In  the  following  areas: 


a.  Is  the  report  adequate  to  support  the  development  of  a  test 
program  to  further  research  pressure  vessel  failures? 

b.  What  Improvements  would  make  the  report  more  adequate? 

c.  Are  there  any  general  comments  concerning  the  report? 


3.  Please  note  any  specific  errors  which  have  been  discovered. 
Include  the  page  number  for  reference . 


V 


PRSTACB 


Work  on  the  study  leading  to  this  report  began  in  order  to  identify 
whether  any  methodologies  exist  for  determining  the  hazards  associated 
with  a  pressure  vessel  failure  and  the  safe  location  of  pressure  vessels 
with  respect  to  nearby  facilities.  As  more  references  were  reviewed,  it 
became  clear  that  most  of  the  test  data  on  blast  waves  and  fragments 
were  specific  to  munitions  blasts  and  did  not  model  the  failure  of  a 
pressure  vessel  very  closely.  Most  of  the  theoretical  work  and  computer 
codes  developed  specifically  for  vessel  bursts  have  not  been 
sufficiently  verified  due  to  a  lack  of  experimental  data.  As  a  result, 
a  major  goal  of  this  study  became  the  identification  of  additional 
testing  necessary  to  validate  methodologies  to  predict  blast  wave 
overpressure  and  fragment  initial  velocity  and  range. 

Program  work  continues  with  the  development  of  a  series  of 
comprehensive  test  plans  for  vessel  bursts,  using  the  Preliminary  Test 
Program  contained  herein  as  a  starting  point.  Initial  conduct  of  vessel 
bursts  is  planned  for  the  summer  of  1989  and  continuing  into  subsequent 
years.  Vessels  and  instrumentation  are  currently  being  sought  for  the 
tests.  Interested  parties  are  invited  to  contact  the  authors. 
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1 . 1  Report  Summation 

The  hazards  associated  with  the  failure  of  pressure  vessels  are  of 
considerable  Interest  to  those  responsible  for  safeguarding  against  such 
failures.  In  general,  much  study  has  been  afforded  toward  the 
identification  of  pressure  vessel  failure  modes  and  has  resulted  in  the 
development  of  methods  for  the  design  against  failure;  such  as  the 
development  of  the  ASME  Boiler  and  Pressure  Vessel  Code,  Section  Vlll. 
In  some  applications,  a  prediction  of  the  hazards  of  a  failed  pressure 
vessel,  and  the  effects  on  its  surroundings,  is  warranted.  Methods  for 
predicting  hazards  of  pressure  vessel  failures,  however,  have  received 
much  less  study  and  a  consensus  methodology  is  not  to  be  found. 
Establishing  the  potential  hazard  is  required  for  the  effective 
mitigation  of  the  consequences  of  failure  through  facility  design  and 
operation . 

The  objective  of  this  study  is  to  identify  the  current  methods 
available  to  establish  the  hazard  a.S3ociated  with  the  failure  of  a 
pneumatic  pressure  vessel  {see  Figure  1-1) .  It  is  desirable  to  safely 
locate  the  vessels  with  respect  to  nearby  facilities,  systems  and 
equipment .  This  involves  determining  the  range  at  which  blast 
overpressure  and  fragments  jeopardize  the  integrity  of  the  facilities, 
systems  and  equipment . 

The  approach  used  in  this  study  was  to  review  the  available 
literature  on  the  theory,  experimental  data  and  computer  codes  for 
determining  the  effects  of  a  blast  wave  and  characteristics  of  the 
fragments  generated  as  a  result  of  the  burst. 

Methodologies  are  considered  to  predict  the  overpressure  from  a 
blast  wave  and  the  initial  velocities  and  ranges  of  fragments  produced 
from  the  vessel  burst.  The  major  limitation  is  that  no  methodology  has 
been  adequately  validated  by  experimental  data. 
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Figure  1-1.  Report  Sutnnation 
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The  report  concludes  with  a  recommendation  for  testing  to  generate 
the  data  necessary  to  verify  methodologies  for  predicting  pressure 
vessel  burst  hazards.  A  preliminary  test  program  is  presented  to 
address  this  concern. 

The  scope  of  this  report  is  limited  to  the  determination  of  the 
blast  wave  and  the  determination  of  fragmentation;  or  specifically,  the 
overpressure  in  air  versus  distance  from  the  gas  filled  pressure  vessel 
rupture,  and  the  quantity  of  vessel  material  pieces  following  rupture 
and  their  respective  velocities  and  masses.  The  studied  rupture  is 
without  detonation,  deflagration  or  other  chemical  reaction.  The 
effects  of  the  generated  blast  wave  and  fragmentation  on  the 
surroundings  (e.g.,  facilities,  personnel)  and  the  generation  of 
secondary  fragmentation  or  debris  is  not  discussed. 

1.2  Current  Blast  Wave  Determination 


Although  limited,  there  have  been  milestone  studies  developed  and 
published  on  the  prediction  of  blast  waves  generated  by  pressure  vessel 
ruptures.  These  were  found,  however,  to  differ  in  methods  and  results 
and  thereby  led  to  the  development  of  this  study  and  report.  These 
studies  are  discussed  in  Section  2  with  specific  emphasis  given  to  those 
factors  that  can  vary  the  determination  of  blast  wave  strength  versus 
distance,  such  as  expansion  energy  available  in  a  compressed  gas,  and 
the  comparison  of  a  vessel  rupture  with  TNT  detonation.  In  addition, 
and  perhaps  more  importantly,  the  results  of  a  literature  search  for 
actual  pressure  vessel  burst  test  data  are  discussed  and  compared  to 

current  blast  wave  prediction  methods.  A  brief  review  of  computer 

models  developed  for  blast  wave  prediction  is  also  presented. 

1.2.1  Bac)cground 

A  pressure  vessel  rupture  results  in  the  uncontrolled  release  of 
the  potential  energy  of  the  compressed  gas  to  the  surroundings.  A 
rupture  of  small  size  (i.e.,  leak),  with  low  tearing  force,  will 

dissipate  the  stored  energy  in  the  vessel  but  over  a  relatively  long 
period  of  time  [1]  and  may  not  generate  a  blast  wave.  An  explosive 
rupture  of  a  pressure  vessel,  where  the  stored  energy  is  released 
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instantaneously/  wtuld  create  a  blast  wave  (i.e.,  shockwave)  in  the 
surrounding  air.  For  the  purpose  of  analysis  of  available  expansion 
energy  and  prediction  of  blast  waves,  the  vessel  material  is  assumed  to 
completely  disintegrate  upon  rupture  with  no  work  performed  on  the 
vessel  fragments.  However,  10%  or  more  of  the  available  energy  is 
expended  to  deform  the  vessel  structure  and  impart  velocity  to 
fragments . 

Explosive  disintegration  will  generate  a  blast  wave  resulting  in 
overpressure  (pressure  above  atmospheric)  at  the  vessel  surface  equal  to 
the  pressure  in  the  vessel  [1] .  As  the  blast  wave  advances,  the  energy 
is  spread  over  the  waves'  spherical  frontal  area,  this  area  increasing 
with  the  square  of  the  distance  from  the  point  of  rupture  [2] . 
Overpressure,  blast  wave  velocity  and  therefore,  blast  effect,  decrease 
rapidly  with  distance.  After  passage  of  the  shockwave,  the  pressure 
decreases  progressively  until  a  suction  phase  follows  in  which  pressure 
drops  below  normal  atmospheric  pressure.  The  negative  pressure  is  a 
result  of  the  spherical  outrush  of  gases  from  the  center  of  the  rupture 
cau.sing  an  overexpansion  (4]  .  The  pressure  above  atmospheric  at  the 
shock  wave  front  is  the  peak  overpressure  and  is  the  blast  wave 
characteristic  most  often  used  to  establish  the  relative  hazard  (i.e., 
shock  wave  intensity)  associated  with  ruptures  and  explosions  at  a  given 
distance.  (The  potential  hazard  to  personnel  and  structures  for  varying 
overpressures  is  very  well  documented.)  Overpressure  is  related  to 
other  blast  wave  characteristics  such  as  velocity  and  impulse,  and  is 
usually  presented  as  ".side-on"  pressure,  that  is,  the  overpressure  which 
would  be  observed  when  there  is  no  interaction  between  blast  and  the 
structure.  A  blast  wave  incident  with  a  structure,  such  a?  a  wall, 
generates  a  reflected  shockwave  This  initial  reflected  pressure  may  be 
several  times  as  great  as  the  side-on  peak  overpressure. 

1.2.2  Peak  Overpressure  Determination 

Determining  the  total  amount  of  energy  transferred  from  a  pressure 
vessel  rupture  (or  any  explosion)  to  the  resulting  blaatwave  is  required 
to  assess  explosive  blast  and  peak  overpressure.  The  explosive  anergy 
from  the  rapid  expansion  of  compressed  gas  can  be  determined  by 
application  of  basic  thermodynamrc  relationships  that  are  a  function  of 
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prassure,  volume,  and  temperature.  The  expansion  is  most  often  assumed 
to  be  tsentropic,  defined  as  adiabatic  (no  heat  transfer)  and  reversible 
(no  losses  due  to  friction,  intermixing  etc.).  However,  isothermal 
expansion  is  also  considered  applicable  by  some  references  (for  example, 
the  Air  Force  System  Command  Design  Handbook)  for  the  determination  of 
expansion  energy.  Isothermal  expansion  implies  heat  is  transferred  into 
the  expanding  gas,  maintaining  constant  temperature.  Assuming  this 
expansion  process  results  in  a  calculated  energy  of  approximately  three 
times  greater  than  for  isentropic  expansion.  Section  2  derives 
relationships  for  both  isentropic  and  isothermal  expansion  energy  for 
ideal  gas  and  discusses  the  applicability  of  each.  The  influence  of 
expanding  real  versus  ideal  gas  is  also  presented. 

By  assuming  the  expansion  energy  of  the  initially  compressed  gas 
develops  a  blastwave  and  that  the  blast  generated  from  a  pressure  vessel 
rupture  is  the  same  as  the  blast  generated  by  a  TNT  (symmetrical  2,  4, 
6-trinitrotoluene)  explosion,  the  peak  overpressure  versus  distance  for 
a  pressure  vessel  rupture  could  then  be  taken  as  that  of  a  TNT  explosion 
of  equivalent  blast  energy  yield.  This  approach  is  frequently  used, 
since  significant  data  exists  on  TNT  explosions  and  the  resulting  blast 
wave . 


The  data  is  usually  presented  as  peak  overpressure  versus  distance 
per  unit  weight  of  TNT.  Discussion  of  converting  expansion  energy 
(e.g.,  ft-lbs)  to  equivalent  weight  of  TNT  is  provided  is  Section  2.3. 

1 . 3  Current  Fragmentation  Determination 

There  has  been  some  study  in  the  last  20  years  on  the  prediction  of 
characteristics  of  fragments  generated  from  pressure  vessel  ruptures. 
These  studies  were  found  to  be  generally  related  to  the  same  theory  for 
establishing  the  energy  imparted  to  the  fragments  and  varied  principally 
in  the  assumptions  necessary  to  model  a  specific  vessel  rupture. 
Further,  as  with  blast  wave  determination,  limited  testing  has  been 
accomplished  toward  predicting  fragmentation  characteristics.  Published 
studies  on  fragments  from  bursting  pressure  vessels  are  discussed  in 
Section  3.0  with  specific  purpose  given  to  those  factors  that  determine 
the  fragment  hazard,  such  as  velocity  and  distance  traveled.  The 


results  of  actual  pressure  vessel  burst  test  fragmentation  data  are  also 
discussed,  as  well  as  a  brief  review  of  computer  codes  developed  from 
fragment  characteristics. 

1.3.1  Background 

A  pressure  vessel  rupture  occurs  with  structural  failure  of  the 
vessel  material.  Failure  is  due  to  an  over stressed  condition  resulting 
in  crack  initiation  and  propagation  caused  by  either  defects,  designed 
stress  concentrations,  missile  impact  or  overloading.  If  crack 
propagation  is  stable  such  that  the  overstress  is  relieved  following 
crack  growth,  a  rupture  of  small  size  (x.e.  leak)  may  occur  and 
dissipate  the  stored  energy  in  the  vessel  without  the  generation  of 
fragments.  When  crack  propagation  is  unstable,  an  explosive  rupture  of 
the  vessel  may  occur  with  the  vessel  material  separating  into  two  or 
more  fragments.  The  expanding  gas  will  accelerate  those  fragments  no 
longer  integral  with  the  vessel  or  otherwise  fixed. 

To  estimate  the  hazard  of  fragments,  regardless  of  source,  complete 
knowledge  of  the  fragment  characteristics  is  required;  this  includes 
fragment  mass,  shape,  drag  coefficient  and  initial  velocity  (speed  and 
direction)  .  For  ground-based  vessels,  this  data  can  be  used  with 
ballistic  type  calculations  to  determine  the  impact  range  and  terminal 
kinetic  energy.  However,  determining  these  fragmentation 
characteristics  is  exceedingly  difficult.  The  most  difficult  problem  is 
the  prediction  of  the  total  number  of  fragments  and  the  mass  of  each 
individual  fragment.  The  current  state-of-the-art  analysis  of 
fragmentation  does  not  provide  any  analytical  theory  to  predict  either 
the  individual  mass  or  the  total  number  of  fragments. 

1.3.2  Fragment  Velocity  and  Range  Determination 

Fragments  produced  upon  vessel  burst  are  accelerated  away  from  the 
vessel  by  the  gas  remaining  within  the  vessel .  The  pressure  of  the 
still  contained  gas  decreases  as  gas  escapes  the  confines  of  the  vessel 
remains.  As  this  gas  pressure  decreases,  the  fragment  acceleration  also 


decreases.  The  gas  pressure  can  be  written  in  terms  of  an  appropriate 
equation  of  state  for  either  an  ideal  or  real  gas.  Initial  fragment 
velocity  is  determined  when  the  acceleration  becomes  negligible. 

Fragment  range  is  a  function  of  the  fragment  trajectory.  The 
trajectory  is  dependent  upon  the  fragment's  initial  speed  and  launch 
angle,  mass  and  coefficient  of  drag  (or  lift)  while  in  flight.  The 
coefficient  of  drag  or  lift  depends  upon  the  speed,  size,  and  shape  of 
the  fragment  and  can  generally  be  found  from  tables  or  estimated. 
Fragment  range  would  then  be  determined  by  a  computer  program  using 
numerical  integration  of  the  basic  ballistic  equations. 

1 . 4  Bibliography 

An  exhaustive  bibliography  of  literature  pertinent  to  blast  waves 
and  fragmentation  is  contained  in  Reference  (6)  .  A  number  of  these 
publications  are  applicable  to  the  failure  of  pneumatic  pressure  vessels 
and  were  reviewed  during  the  preparation  of  this  report.  Those  which 
were  finally  used  as  sources  for  this  report  are  listed  in  the 
References  herein. 

1 . 5  Summary  of  Conclusions  and  Recommendations 

The  two  major  methodologies  for  predicting  blast  wave  overpressure 
involve  isothermal  and  isentropic  expansion  of  a  gas  which  might  be 
either  real  or  ideal.  Both  models  tend  to  be  overconservative  in  the 
near  field  in  that  they  predict  greater  overpressures  than  are  actually 
observed  Of  tne  two  methodologies,  the  isentropic  model  appears  to  be 
the  more  accurate  and  suitable  for  verification.  The  ideal  gas 
approximation  is  also  overconservative;  either  the  ideal  or  real  gas 
model  might  be  suitadjle  for  further  study.  A  lac)c  5f  available  blast 
data  specific  to  pressure  vessels  makes  impossible  the  verification  of  a 
methodology  of  either  of  two  comput«,r  codes  currently  available. 

There  are  several  methodologies  available  to  predict  initial 
fragment  velocity;  these  models  differ  principally  in  their  as.iumptions 
and  applications.  Two  limitations  generally  apply  to  all  of  these 
models:  (1)  input  data  is  required  that  involves  knowledge  of  how  the 
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vessel  fails  or,  at  least,  what  fragments  are  generated,  and  (2)  there 
is  insufficient  data  to  verify  any  of  the  methodologies.  The  available 
computer  codes  suffer  from  the  same  limitations. 

It  is  much  easier  to  predict  the  ranges  of  the  fragments,  given 
their  initial  velocities.  Computer  codes  are  also  available. 

It  is  strongly  recommended  that  additional  testing  be  conducted  to 
obtain  the  necessary  data  to  validate  the  blast  wave  and  fragment 
methodologies.  A  preliminary  test  program  is  presented  to  provide 
guidance  in  setting  up  testing. 
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SECTION  2.  BLAST  NAVE  THEORY  AND  EXPERIMENTAL  APPROACHES 


The  purpose  of  this  section  is  to  describe  current  methodologies 
that  analyze  the  effects  of  a  blast  wave  generated  by  the  burst  of  a 
pneumatic  vessel.  The  section  begins  with  an  introductory  description 
of  the  blast  wave  formed  by  a  vessel  burst.  Four  models  are  then 
developed  to  approximate  the  energy  release  in  the  blast  wave;  an  ideal 
gas  is  assumed  to  expand  isentropically,  adiabatically  and  isothermally . 
The  isentropic  and  isothermal  models  are  chosen  for  further  review  and 
compared.  Next,  the  differences  between  real  and  ideal  gas  are 
discussed  and  the  isentropic  and  isothermal  models  are  compared  using 
ideal  gases.  The  energy  release  from  a  bursting  vessel  is  then  equated 
to  a  TNT  explosion,  and  the  TNT  equivalency  is  compared  with 
experimental  data.  Finally,  computer  codes  which  have  been  gsnerated  to 
predict  the  blast  waves  from  vessel  bursts  are  reviewed. 

2 . 1  Theory 

The  bursting  of  a  high  pressure  vessel  results  in  an  instantaneous 
release  of  the  energy  stored  within  the  vessel.  The  energy  is 
dissipated  through  various  means,  including  rupture  of  the  vessel, 
propulsion  of  fragments,  and  a  blast  wave. 

In  general,  for  failures  of  vessels  filled  with  a  compressed  gas, 
it  may  be  assumed  that  the  gas  expands  in  an  adiabatic  manner.  Op  to  an 
initial  pressure  of  about  IbOU  psi  (about  100  atmospheres)  it  can 
further  be  assumed  that  the  gas  expands  in  an  ideal  manner.  Above  1500 
psi,  the  ideal  gas  laws  are  not  accurate  and  real  gas  effects  become 
important,  causing  the  actual  energy  release  to  be  less  than  that 
predicted  by  ideal  gas  laws  f3] .  Therefore,  the  use  of  ideal,  adiabatic 
expansion  equations  for  vessel  pre.ssures  above  about  100  atmospheres 
will  produce  conservative  results.  Ideal  gas  and  real  gas  energy 
releases  are  examined  later  in  this  section. 

The  mediums  in  pressure  vessels  can  be  put  into  three  categories 
from  the  highest  hazard  potential  from  a  blast  wave  to  the  least: 

(1)  compressed  gases,  (2)  compressed  flash-evaporating  liquids,  and  (3) 
compressed  liquids.  These  hazard  categories  are  based  on  the  stored 
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energy  per  unit  volume.  Pressure  vessel  failure  can  occur  for  numerous 
reasons  with  the  major  ones  being  the  following:  overpressurization, 
material  failure,  and  missile  impact. 

This  report  analyzes  the  energy  stored  in  a  compressed  (ideal  and 
real)  gas  pressure  vessel  based  upon  the  thermodynamic  relationships  for 
isentropic,  adiabatic,  and  isothermal  expansions. 

2.1.1  Introduction  to  Blast  Wave  Characteristics 

A  blast  wave  (shock)  is  formed  when  the  atmosphere  surrounding  a 
vessel  that  fails  is  forced  away  by  the  expanding  fluid  from  the  vessel . 
Figures  2-la  &  lb  illustrate  theoretical  pressure  vessel  blast  waves  in 
the  form  of  a  pressure-time  curve  at  a  fixed  location  from  the  vessel. 
Some  experts  say  the  blast  wave  is  characterized  by  a  sharp  rise  to  a 
peak  pressure  and  an  exponential  decay  thereafter  (Figure  2-la)  while 
others  contend  that  there  is  a  less  rapid  rise  to  a  less  distinct  peak 
with  a  gradual  decay  to  ambient  pressure  (Figure  2-lb) .  As  more 
experimental  pressure  vessel  burst  data  is  generated,  a  better 
determination  as  to  blast  wave  characteristics  will  be  obtainable. 
Figure  2-2  shows  a  theoretical  pressure-distance  plot  for  a  blast  wave 
at  .several  specific  times.  Once  any  blast  wave  has  traveled  several 
vessel  diameters  from  the  source,  it  tends  to  follow  these 
configurations.  Near  the  vessel,  there  can  be  differences  in  the  curves 
depending  on  the  source  of  the  blast . 

It  can  be  seen  from  Figure  2-la  that  there  are  both  positive  and 
negative  pressure  phases.  In  general,  the  positive  pressure  pulse  is 
much  stronger  and  has  greater  effects  than  the  negative  pulse  [4]  .  An 
important  point  to  notice  from  Figure  2-2  is  that  as  the  distance  from 
the  blast  wave  source  increases,  the  peak  overpressure  decreases. 
Eventually,  the  peak  overpressure  will  decrease  to  a  point  where  no 
damage  will  be  done  to  structures,  equipment,  or  personnel.  Methods 
will  be  described  later  to  relate  distance  to  overpressure. 


10 


PRESSURE 


Overpressure  is  an  important  factor  in  determining  safe  separatic  ■ 
distances,  because  once  the  overpressure  vs.  distance  is  determined,  it 
is  possible  to  enter  a  table  that  identifies  potential  damage  from 
overpressures.  (Note  that  there  are  various  such  tables  published,  many 
of  which  give  different  data.) 

Referring  bac)c  to  Figure  2-la,  it  can  be  seen  that  ac  a  specific 
distance  from  the  blast,  the  blast  wave  arrives  as  an  Instantaneous 
pressure  increase,  which  immediately  begins  to  decay  and  eventually 
becomes  negative.  Prior  to  shock  front  arrival,  the  pressure  is  ambient 
pressure  P^.  At  arrival  time  t.,  the  pressure  rises  abruptly  to  a  pea)c 
value  P,*+  P,.  The  pressure  then  decays  to  ^UQble^t  in  total  time  t,  +  T*, 
drops  to  a  partial  vacuum  of  amplitude  P.',  and  eventually  returns  to  P„ 
in  total  time  t,  +  T*  +  T‘  The  quantity  P,*  is  usually  termed  the  pea)c 
side-on  overpressure,  or  merely  the  peak  overpressure.  The  portion  of 
the  time  history  above  initial  ambient  pressure  is  called  the  positive 
phase,  of  duration  T* .  That  portion  below  P„,  of  amplitude  P,'  and 
duration  T'  is  called  the  negative  phase  (5) . 

Now  that  these  basic  concepts  of  a  blast  wave  have  been  introduced, 
the  theoretical  approaches  to  determining  the  energy  released  by  the 
failure  of  a  gas-filled  pressu- e  vessel  will  be  discussed.  A  critical 
problem  has  been  to  accurately  assess  the  energy  release  as  a  result  of 
the  accident  or  incident.  Current  prediction  methods  determine  the 
energy  stored  in  the  pressure  vessel  and  equate  it  Vo  an  equivalent 
quantity  of  TNT  (equivalency  based  on  energy  release) .  This  method  is 
used  because  a  large  amount  of  experimental  data  exists  for  blast  waves 
generated  from  TNT  explosions.  Although  the  comparison  with  TNT  is 
convenient,  it  will  be  shown  in  Section  2.3  that  the  correlation  is  at 
best  a  fair  approximation. 

2.1.2  Energy  Released  by  Isentropic  Expansion  of  an  ideal  Gas 

The  blast  eneigy  is  defined  as  the  energy  released  by  the 
pressurized  gas  following  vessel  failure.  In  this  subsection,  a 
relationship  for  the  blast  energy  based  upon  isentropic  expansion  of  an 
ideal  gas  will  be  derived,  beginning  with  the  fir.*it  law  ot 
thermodynamics . 
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The  system  of  concern  in  this  study  consists  of  gas  under  high 
pressure  contained  in  a  pressure  vessel  at  a  ground  level .  Initially, 
this  system  is  at  pressure  P,  and  it  contains  a  volume  V,  of  gas.  After 
vessel  rupture,  the  gas  in  the  vessel  expands  to  P,  -  P.^j,  and  to  an 
unknown  volume  v, . 

The  first  law  of  thermodynamics  can  be  written  as: 

q-w-Ae-Au  +  Ake  +  Ape  ■••...  (i) 

where  :  Q  -  heat  transferrea 

W  ••  work  done  on  surroundings 
AE  -  change  in  total  energy 
Au  -  change  in  internal  energy 
Ake  ■  change  in  kinetic  e.nergy 
Ape  -  change  in  potential  energy 

...  -  change  in  other  energies  (e.g.,  chemical,  electrical) 

Because  the  expansion  is  isentropic,  (i.e.,  adiabatic  and 
reversible),  Q  -  0  (Later  sections  will  examine  other  processes.). 
Because  the  system  is  initially  at  rest,  and  is  also  at  rest  following 
the  expansion  process,  Ake  is  negligible.  Similarly,  Ape  is  negligible 
because  we  shall  assume  the  entire  process  takes  place  at  ground  level. 
Other  energies  such  as  chemical  and  electrical  are  not  applicable  to 
this  process  and  are  also  negligible.  These  assumptions  reduce  Equation 
(1)  to: 


-W  -  Ae  -  Au 

The  ideal  gas  law  states  that  for  the  expansion  of  a  gas: 
W  -  -C,  At 

where:  ••  constant  volume  specific  heat 

At  -  Tj  -  T;  (temperature  change) 


(2) 


(3) 
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The  constant  volume  specific  heat  is  a  property  of  a  fluid  and  can  be 
written  as: 


C.  - 

K-1  (4) 


where:  R  -  universal  gas  constant 

K  -  specific  heat  ratio  (for  a  particular  ideal  gas) 


Inserting  Equation  (4)  into  Eq^iation  (3),  the  relationship  for  v/or)t 
becomes : 

W  -  _R_  (T,  -  T,) 

K-1  (5) 


The  basic  ideal  gas  law  is: 

PV  -  RT 
or 

R  -  p_y 

T  ' 

Substituting  Equation  (7)  into  Equation  (5)  and  reducing: 

W  -  (1  -  T,) 

K-1  T,  (8) 


(6) 

(7) 


An  isentropic  ideal  gas  relation  is: 


K-1 

I;  -  1:  K  (9) 

T,  P; 


Substituting  Equation  (9)  into  Equation  (8): 


W  - 


K-1 

K 


(10) 
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Equation  (10)  gives  the  isenttopic  energy  released  by  the  failure 
of  a  vessel  containing  a  volume  of  ideal  gas,  V,,  at  a  pressure  of  P, . 
P,  is  the  surrounding  atmospheric  pressure.  From  this  equation,  the 
energy  available  from  a  given  vessel  can  be  calculated.  The  isentropic 
formula  is  used  often  in  gas  dynamics  because  it  is  a  simple  formula, 
and  the  friction  losses  and  dissipation  losses  are  usually  very  minimal. 

The  specific  heat  ratio  K  is  essentially  independent  of  pressure 
and  temperature  for  monatomic  gases.  For  diatomic  gases  and  other 
gases,  specific  heat  ratio  decreases  somewhat  with  increasing 
temperature.  Data  are  tabulated  and  readily  available  (24)  if  increased 
accuracy  is  desired. 

As  will  be  seen  later  in  Section  2.3,  this  energy  can  be  converted 
into  an  equivalent  quantity  of  liiT.  Then,  the  blast  wave  associated 
with  the  vessel  failure  can  be  assumed  to  be  equivalent  to  the  effects 
of  that  quantity  of  TNT. 

In  a  real  situation,  the  blast  wave  energy  released  by  the  vessel 
failure  would  be  less  than  that  of  the  theoretical  value  obtained 
through  Equation  (10)  .  One  reason  for  this  has  been  discussed  -  at 
pressures  above  about  1500  psi  and  ambient  temperature,  real  gas  effects 
become  significant  and  reduce  the  energy  released  as  compared  to  an 
ideal  gas.  (Real  gas  effects  are  discussed  later  in  this  section.) 
Another  reason  that  Equation  (10)  or  the  equations  developed  in  the  next 
two  subsections  would  give  a  conservative  result  is  that  a  portion  of 
the  energy  contained  in  the  piessu:  zed  gas  would  not  go  into  the  blast 
wave,  but  instead  would  oe  axp  ded  through  vessel  rupture  and 
fragmentation.  Estimates  of  the  energy  expended  through  a  shoOc  wave 
range  from  20%  to  80%  of  the  total  energy  available  [6]  .  In  real 
situations,  the  energy  expended  through  vessel  rupture  and  fragmentation 
will  vary  based  upon  material  characteristics  (e.g.,  fracture  toughness, 
material  flaws,  and  mode  of  vessel  failure) . 
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2.1.3  Energy  Released  by  Adiabatic  Expansion  of  an  Ideal  Gas 

This  subsection  develops  the  energy  release  equation  based  upon 
adiabatic  expansion,  whereas  the  previous  subsection  examined  isentropic 
(i.e.,  adiabatic  and  reversible)  expansion. 

Again  starting  with  the  first  law  of  thermodynamics: 

Q-w  -  Ae  “  Au  +  Are  +  Ape  +  . . .  (ii) 

As  was  the  case  for  the  previous  section,  Q  “  0,  and  Ake  and  Ape  are 
negligible. 

The  development  of  the  adiabatic  equation  is  identical  to  the 
development  of  the  isentropic  equation  up  through  Equation  (8) : 


At  this  point,  the  temperatures  can  be  replaced  with  the  following 
relationship : 


T  -  PV  (13) 

R 


Substituting  Equation  (13)  into  Equation  (12),  and  reducing: 


”  “  LYi  /  ^  I 

K-1  p.v.y 


(14) 


Also,  the  total  volume,  V,  equals  the  specific  volume,  v,  multiplied  by 
the  mass,  m; 


V  -  mv 


(15) 


Also,  V  - 

P 


where  p 


density 


(16) 
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Substituting  Equations  (15)  and  (16),  Equation  (14)  now  becomes: 


Equation  (17)  gives  the  maximum  wor)c  that  gas  undergoing  an  adiabatic 
expansion  can  perform  on  its  surroundings. 


Adiabatic  means  that  no  heat  is  transferred  into  or  out  of  the 
system.  This  is  the  most  likely  situation  during  tank  burst,  because 
the  gas  expands  so  fast  it  doe.s  not  have  time  to  transfer  significant 
heat  and  the  adiabatic  formula  does  not  assume  the  expansion  to  be 
reversible.  However,  this  formula  is  not  widely  used  because  it 
requires  knowledge  of  the  gas  density  at  point  2. 


2.1.4  Energy  Released  by  Isothermal  Expansion  of  an  Ideal  Gas 

This  subsection  develops  the  energy  release  equation  based  upon 
isothermal  expansion,  whereas  the  previous  two  subsections  examined 
isentropic  and  adiabatic  expansion.  The  isothermal  process  can  be 
modeled  by  adding  a  heat  source  to  the  isentropic  process  of  subsection 
2.1.2,  so  that  heat  would  be  available  to  the  system  during  the 
expansion  process  making  temperature  constant. 

For  an  isothermal  process,  T,  -  T,.  Starting  with  the  first  law  of 
thermodynamics : 

Q  -  w  -  Ae  -  Au  -i  Ake  +  Ape  (18) 


In  this  case  Q  ¥  0  because  of  the  heat  source.  Again,  Ake,  Ape  and 
changes  in  other  energies  are  negligible,  so  Equation  (18)  reduces  to: 

Q  -  W  -  Au  (19) 

But  Au  is  defined  as: 
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Because  the  process  is  isothermal,  dT  -  0;  therefore,  Au  -  O,  so 
Equation  (19)  is  further  reduced  to: 


Q-W-OorW-Q 


(21) 


Wor)c  can  be  defined  in  terms  of  pressure  and  volume  as  follows: 

«  -/pdv  (22) 


And  since  PV  -  RT,  or  P  ■  ^  : 

V 


W  -  /*RT 

J  V 


dv 


(23) 


R  is  the  universal  gas  constant.  Because  we  are  examining  an 
isothermal  process,  T  is  also  constant. 


Equation  (23)  now  becomes: 


W  «  RT 


W  -  RT  In  Vj 

V, 


(24) 


as : 


The  ideal  gas  law  of  PV  -  RT  can  be  written  for  two  separate  points 


P,V,  P,V,  and  since  Tj  -  T, 

T.  T, 


P,V,  =  P,V,,  or 


P,  -  V, 
P.  V, 


(25) 


Substituting  Equation  (25)  into  Equation  (24): 


W  -  RT  In  P, 
P. 


(26) 
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Also,  PiV,  -  RT 


So  Equation  (26)  now  becomes: 

W  -  P,Vi  In  P^  (27) 

P: 

P,,  V,,  and  P,  are  all  determinable;  therefore,  the  maximum  wor)^ 
that  an  isothermal  expansion  can  do  to  its  surroundings  can  be  calcu¬ 
lated  using  Equation  (27) . 

2.1.5  Comparison  of  Energy  Releases  Calcu'ated  by  Isentropic,  Adiabatic 
and  Isothermal  Expansion  Equation^  Ideal  Gas 

Equations  (10),  (17),  and  (27)  present  the  isentropic,  adiabatic 

and  isothermal  expansion  formulas,  respectively,  for  ideal  gas.  This 
section  will  compare  the  results  of  the  isentropic  and  isothermal 
equations  for  vessel  failure  where  initial  pressures  (P,)  range  from  100 
to  15,000  psi  and  Pj  is  assumed  to  be  atmospheric  pressure  (14.7  psia)  in 
all  cases. 

The  isothermal  expansion  formula  is  a  single  parameter  (pressure) 
formula,  easy  to  use.  However,  for  the  non-burning  high  pressure  burst, 
isothermal  expansion  is  not  realistic,  due  to  the  fact  that  there  is  no 
heat  generated  to  support  constant  temperature  gas  expansion.  For  this 
reason  TNT  equivalency  based  on  isothermal  expansion  is  over¬ 
conservative.  The  adiabatic  expansion  formula  requires  both  tan)^ 
pressure  and  density  to  determine  the  TNT  equivalency.  Adiabatic  means 
no  hear,  transfer  in  or  out  of  the  system.  This  is  the  most  li)tely 
situation  during  tank  bursts,  because  the  gas  expands  so  fast  it  does 
not  have  time  to  transfer  significant  heat  into  or  out,  and  the 
adiabatic  case  does  not  assume  the  expansion  to  be  reversible.  This 
formula  is  not  widely  used  because  it  requires  knowledge  of  the  gas 
density  at  point  2.  The  isentropic  relationship  reduces  the  two- 
parameter  adiabatic  formvila  to  a  one-parameter  (pressure)  formula. 
However,  there  is  concern  over  the  fact  that  isentropic  implies 
reversible  because  the  reversible  process  is  an  idealization  [7] .  It  is 
a  concept  which  can  be  approximated  very  closely  at  times,  but  never 


matched.  In  many  actual  cases,  effects  such  as  friction,  electrical 
resistance,  and  inelasticity  can  be  substantially  reduced,  but  their 
complete  elimination  is  not  found.  These  and  similar  effects  are 
frequently  called  dissipative  effects,  since  in  all  cases  a  portion  of 
the  energy  originally  in  the  system  is  converted  or  dissipated  into  a 
less  useful  form.  Only  in  the  absence  of  dissipativ  effects  can 
certain  forms  of  energy  be  converted  into  other  forms  without  any 
apparent  loss  in  the  capabilities  of  the  system.  Any  system  which  is 
returned  to  its  initial  state  after  experiencing  an  irreversible  process 
will  leave  a  history  in  the  surroundings  due  to  irreversibilities.  A 
partial  list  of  these  effects  relating  to  pressure  vessel  bursts  are: 
shock  waves,  unrestrained  expansion  of  a  fluid,  mixing  of  dissimilar 
gases,  and  mixing  of  the  same  two  fluids  initially  at  different 
pressures  and  ten^eratures .  Nevertheless,  the  reversible  processes  have 
been  found  to  be  an  appropriate  starting  place  on  which  to  base 
engineering  calculations.  This  is  why  the  isentropic  formula  is  the 
most  often  used  formula  in  gas  dynamics. 

Table  2-1  shows  calculated  "alues  for  equations  (10)  and  (27)  . 
Figure  2-3  shows  plots  of  energy  versus  vessel  pressure  for  a  comparison 
of  isentropic  and  isothermal  expansion  formulas.  Energy  was  calculated 
per  unit  volume  for  purposes  of  comparison. 

2 . 2  Real  Gas  Theory 

The  ideal  gas  assumption  for  high  pressure  (greater  than  15C0  psi) 
ruptures  gives  expansion  energies  that  are  unrealistically  high. 
Accurate  estimates  of  blast  parameters  from  high  pressure  bursts  require 
calculations  based  on  real  gas  equations  of  state  supported  by  empirical 
data.  Ideal  gas  behavior  is  adequate  for  most  low  pressure  situations 
(<1500  psi) .  The  expansion  energy  from  pneumatic  pressure  vessel 
rupture  depends  on  fill  gas  temperature,  ratio  of  specific  heats, 
pressure,  and  volume.  Pittman  [8)  found  in  his  studies  that  for 
monatomic  fill  gas  argon  at  ambient  temperature,  about  77  *F,  increasing 
burst  pressure  above  15,000  psi  results  in  less  significant  increases  in 
expansion  energy  than  at  lower  pressures. 
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Table  2-1 

Comparison  of  Isentropic  and  Isothermal 
Expansion  Energies  for  Ideal  Gas  (GN2) 


Initial  Vessel 

Pressure,  Pj  (psia) 

Isentropic 

Expansion 

Energy,  Eq  (10) 
ft-lb, 
ft' 

isothermal 

Expansion 

Energy,  Eq  (27) 
ft-lb, 
ft' 

0 

0 

0 

50 

8040 

13807 

100 

18334 

33933 

300 

66082 

138840 

500 

118204 

263540 

750 

186278 

435142 

1000 

256347 

618732 

1500 

400235 

1010999 

2500 

696967 

1862027 

5000 

1464214 

4211588 

7500 

2249958 

6748912 

10000 

3046281 

9408315 

2). 


ENERGY 


Figure  2-3.  Comparison  of  Isentropic  and  Isothermal  Expansion  Energies  for  Ideal  GN2 


22 


T^e  effect  of  burst  pressure  on  expansion  energy  is  illustrated  in 
Figure  2-4,  a  graph  of  pressure  versus  specific  volume.  A  path  is  shown 
for  argon,  whereby  it  is  isothermally  compressed,  beginning  at  1 
atmosphere  and  290  K,  to  15,000  psia;  following  vessel  burst,  it 
isentropically  expands.  The  expansion  energy  is  the  work  that  is 
enclosed  on  the  P-v  curve  by  the  path  followed.  If  the  argon  had  been 
isothermally  compressed  to  50,000  psia,  however,  the  additional  decrease 
in  volume  is  very  small.  When  the  gas  isentropically  expands,  the 
increase  in  work  (enclosed  area)  is  correspondingly  small, 

Pittman  went  on  to  determine  that  by  increasing  the  gas 
temperature,  for  a  given  pressure,  the  expansion  energy  will  increase 
until  the  limit  is  reached  where  the  gas  approaches  ideal  behavior,  see 
Figure  2-5.  For  an  ideal  gas,  expansion  energy  is  independent  of 
temperature  and  molecular  weight;  and  the  ideal  gas  internal  energy 
represents  an  upper  limit  to  the  blast  energy  available  from  an  argon 
vessel  rupture.  Pittman  determined  this  occurred  for  argon  at  a 
pressure  of  15,000  psi  as  the  temperature  approached  1970  F.  Higher  gas 
pressures  require  increasingly  higher  temperatures  to  approach  ideal 
behavior.  For  this  study,  vessel  temperature  will  be  ambient  (70  F)  and 
vessel  pressures  wil]  range  from  50  to  15,000  psi.  So  for  vessel 
pressures  <1500  psi  (100  atm)  ,  an  ideal  gas  approximation  will  be  used 
and  for  vessel  pressures  >1500  psi,  a  real  gas  approximation  will  be 
used.  However,  a  problem  arises  for  vessel  pressures  >10,000  psi  in 
that  compressibility  (Z)  data  is  very  limited  for  GAr  and  GN2  in  this 
region.  Yet,  for  GAr  at  ambient  temperature,  increasing  pressure  above 
15,000  psi  only  slightly  increases  expansion  energy.  So  the  problem  of 
limited  (Z)  data  above  10,000  psi  may  not  be  a  significant  problem. 
Further  vessel  burst  testing  may  resolve  this. 

For  real  gases,  especially  at  high  pressures  and  very  low  tempera¬ 
tures,  the  ideal  gas  equation  of  state  (PV  -  RT)  is  not  a  good 
approximation.  There  are  several  equations  of  state  that  can  be  used 
for  real  gases  (e.g.,  Van  der  Waal's,  Beattie-Bridgeman,  Nobel-Abel  and 
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Figure  24.  Argon  P-V  Compression  and  Expansion  Paths  for  Cold  Tanks 


Redlich-Kwong)  using  enqplrical  coefficients  for  the  appropriate  gas 
(25]  .  O.ie  simple  real  gas  equation  of  state  involves  a  compressibility 
factor  Z,  which  is  defined  as: 


Z  “  RT 


Since  RT/P  is  the  id-~:;l  gas  specific  volume  (v  ideal),  the  compres- 
siDility  factor  may  be  considered  a  measure  of  the  ratio  of  the  actual 
specific  volume  to  ideal  gas  specific  volume.  For  an  idea',  gas,  the 
compressibility  factor  is  unity,  but  for  actual  gases  it  can  be  either 
less  than  or  greater  than  unity.  Hence,  the  compressibility  factor 
measures  the  percent  deviation  of  an  ^ctual  gas  from  idea.^  gas  behavior. 
\pplying  this  concept  to  equations  10,  17,  and  27,  we  get  expansion 
energy  relationships  for  real  gaser  ap  follows: 


K-i" 

Isentropic: 

W  -  P, 

V  z(k^) 

..(!;)  ' 

Adiabatic: 

!!  "  L_  i 

V  Z(K-l)  1 

l-P.  P.j 

I 


Isothermal.  W  «  P,  In  Pj 

V  z'  P, 


(28) 


(29) 


(30) 


Eqpiat  ions  28  and  30  present  the  isentrop  .c  ana  Isothermal  expansion 
energies,  respectively,  for  a  real  gas.  This  section  will  compare  the 
results  of  these  equations  for  vessel  failurts  with  initial  pressures 
ranging  from  50  psi  to  high  pressures  typically  »  sed  at  various  facili¬ 
ties.  Table  2-2  shows  calculated  values  for  equations  28  and  30.  The 
lollowing  are  assumptions  made  to  support  the  calculations: 
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1.  Vessels  are  at  ambient  temperature  (540  R) 

2.  P,  is  atmospheric  pressure  (14.7  psia)  in  all  cases 

3 .  Gas :  (^2 

4  .  K  -  1 . 4  for  GN2 

5.  For  6N2:  -  227  R  and  P<,  -  33.5  atm 

6.  P*  ■■  P/Po  and  ••  T/T„.  Absolute  prtssuxes,  and  temperatures 
in  degrees  Kelvin  or  Rankine,  must  be  used. 

The  curves  are  graphed  in  Figure  2-6.  For  purposes  of  comparison,  the 
values  were  calculated  per  unit  volume. 


Table  2-2.  Comparison  of  iseniropic  and  Isothermal  Expansion 
Energies  for  Real  Gas  (GN2) 


Vessel 

Pressure 

(psia) 

z 

Isentropic 

Energy,  Eq  (20) 

ft-lb, 

ft" 

Isothermal 

Energy,  Eq  (30) 

ft-lb, 

~ft^ 

0 

1.0 

0 

0 

50 

1.0 

8040 

13807 

100 

1.0 

18334 

33933 

500 

1.0 

118204 

263540 

750 

. 

188159 

439537 

1000 

.98 

261578 

631360 

1500 

,  97 

412613 

1042267 

2500 

.99 

704007 

1880835 

5000 

1.16 

1262253 

3630679 

7500 

1.37 

1642305 

4926213 

10000 

1.56 

1952744 

6030971 

2 . 3  TNT  Equivalence 

In  Sections  2.1  and  2.2,  the  energy  release  associated  with  a 
pressure  vessel  burst  is  theoretically  predicted  for  the  following 
processes:  isothermal,  adiabatic,  and  isentropic  for  both  ideal  anrl  real 
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g».a9a .  In  this  section,  it  will  be  shown  how  TKT  equivalency  can  be 
used  to  approximate  a  relative  measure  (e.g.,  pounds  of  TNT)  of 
explosive  magnitude  -or  a  bursting  pressure  vessel.  Once  the  TNT 
equivalence  is  determined,  then  the  effects  of  a  pressure  vessel  failure 
can  be  converted  to  overpressure  at  various  distances  from  the  source. 

2.3.1  Determining  TNT  Equivalence 

The  magnitude  of  an  explosion  is  established  by  the  amount  of 
energy  released.  This  can  be  expressed  directly  in  energy  units  such  as 
ft-lb,  joules,  or  calories.  To  express  explosions  numerically  a 
standard,  TNT  (symmetrical  2,4,6  -  trinitrotoluene),  was  established. 
This  is  because  TNT  is  a  chemically  pure  material,  it  is  readily 
available  for  calibration  purposes,  it  is  relatively  safe  to  handle,  and 
for  specimens  of  )tnown  density  and  crystalline  nature,  its  explosive 
effects  are  quite  reproducible  (4).  In  Kinney  [4],  the  standard  pound 
(gram)  of  TNT  is  defined  as  the  energy  of  1.545  x  10*  ft-lb  (4610  J)  . 
This  will  be  the  standard  for  equating  energy  to  equivalent  pounds  of 
TNT  throughout  this  report.  A  further  discussion  of  conversion  factors 
will  appear  in  Section  2.3.2.  The  energy  expansion  equations  derived  in 
Sections  2.1  and  2.2  will  show  how  this  energy  from  a  bursting  pressure 
vessel  is  equated  to  a  TNT  equivalence. 

Example 

Given:  cylindrical  vessel 

hemispherical  heads 
horizontal  geometry 
gas :  GN2 

V,  -  10  ft^  T.  -  227  R  -  33.5  atm 

P,  -  2000  psi 

P,  ^  atmospheric  ••  14.7  psia 
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Assumptions : 


K  -  1.4 

isentropic  expansion 
real  gas 
Find:  TNT  equivalence 


W  -  5.43  X  10*  £t-lb. 


TNT  eq’jivalence  -  3.52  pounds  TNT 

Therefore,  the  energy  from  a  burst  pressure  vessel  of  the  dimen¬ 
sions  given  in  the  example  can  be  equated  to  the  energy  from  a  charge  of 
3.52  pounds  of  TNT.  This  can  be  repeated  for  all  sizes  and  pressures  of 
pressure  vessels.  Figure  2-7  is  a  plot  of  the  relationship  between  TNT 
equivalence  in  pounds  TNT  per  cubic  foot  and  vessel  pressure  in  psi. 
The  plot  shows  two  curves  relating  to  the  energy  expansion  equations 
(e.g.,  isothermal  and  isentropic).  However,  of  concern  at  this  time  is 
the  fact  TNT  energy  may  differ  from  pressure  vessel  failure  energy 
because  not  all  of  the  vessel  energy  is  directed  into  the  blast  wave  as 
is  the  case  for  TNT  energy.  A  significant  portion  of  the  vessel  energy 
may  be  used  to  tear  the  vessel  metal  and  propel  fragments.  This  can 
only  be  resolved  by  further  experimental  work. 

2.3.2  Overpressure  vs.  Distance  for  TNT 

If  the  vessel  burst  energy  is  to  be  equated  to  a  charge  of  TNT,  a 
means  to  convert  pounds  of  TNT  to  pea)t  overpressure  is  needed. 
Overpressure  is  defined  as  the  peak  pressure  above  ambient  conditions. 
Kingery  [9]  contains  a  compilation  of  experimental  blast  parameters 
ver.sus  distance  ‘-.hat  were  measured  on  5  -,  20  100  -,  and  500-ton  TNT 
surface  bursts.  The  data  from  all  four  tests  were  first  processed  to 
obtain  the  "as  read"  values  of  peak  overpressure,  arrival  time,  positive 
duration,  and  positive  impulse.  Cube  root  scaling  and  altitude 
correction.^)  were  applied  to  these  values  to  bring  them  tc-  standard  .sea- 
level  conditions  and  the  equivalent  of  a  one-pound  TNT  charge.  The  peak 
overpressures  were  obtained  from  direct  measurements  provided  by 
pressure  transducers.  The  scale  values  were  then  used  to  determine  the 


30 


POUNDS  TNT/CUBIC  FOOT 


curve  in  Figure  2-8  and  data  in  Table  A-I  of  Appendix  B;  X  -  scaled  dis¬ 
tance  D/W  where  D  -  distance  from  ground  zero  (ft)  and  W  *•  cube 
root  of  TNT  equivalence  (lb*'’)  and  Ap,  -  pea)c  overpressure  (psi)  .  From 
the  TNT  equivalence,  one  can  determine  the  peak  overpressure  at  various 
distances  from  ground  zero  (point  of  vessel  burst) . 

Cube  root  scaling  can  also  be  applied  to  bursts  which  occur  above 
ground  in  order  to,  for  excunple,  minimize  ground  reflection  (26) .  The 
scaling  equation  is  analogous  to  that  used  for  the  horizontal  range. 

X,  -  H/W 

where  ?_  -  scaled  burst  height 
H  -  height  of  burst  (ft) 

W  -  TNT  equivalence  (lbs) 

As  mentioned  earlier  in  Section  2.3.1,  the  standard  for  equating 
vessel  burst  energy  to  equivalent  weight  TNT  will  be: 

1.545  X  10*  ft-lb  energy  -■  one  pound  of  TNT  for  tni.s  report.  This 
value  comes  from  Kinney  (4)  who  gives  the  basis  for  this  value  as 
follows : 

"The  standard  gram  TNT  is  defined  as  the  blast  energy  of 

4610  J  (1100  cal) .  This  dates  to  early  days  of  nuclear 

devices  when  a  standard  ton  TNT  was  defined  as  the  energy 

release  of  one  million  kilocalories." 

Although  we  choo.se  to  use  Kinney's  conversion  factor,  there  is  no 
agreement  in  the  industry  as  to  a  standard.  The  following  are  other 
conversion  factors  used  in  the  industry: 

•  Brown,  S.,  et  al,  'Energy  Release  Protection  -  ASME  - 

SC6000, " 

l:raft  No.  1,  August  1  983. 

1.426  10*  ft-lb  -  1  pound  TMT 
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•  Brown,  S.,  et  al,  "Protection  Against  Pressure  Vessel  And 
Piping  Explosions/Ruptures,"  Tadjle  2. 

(4520  J/g  TNT) 

1.515  X  10‘  ft-lb  -  1  pound  TNT 

•  Pittman,  J.F.,  "Blast  And  Fragment  Hazards  From  Bursting  High 
Pressure  Tanks,"  17  May  1972. 

(1018  Cal/g  TNT) 

1.428  X  10‘  ft-lb  -  1  pound  TNT 

•  Brown,  S.,  "Pressurized  Systems  Energy  Release  Protection," 
NASl-17278,  May  1986. 

(1832.4  Btu/lb  TNT) 

1.425  X  10‘  ft-lb  -  1  pound  TNT 

•  Baker,  W.E.,  et  al,  "Wor)cbook  For  Predicting  Pressure  Wave 
and  Fragment  Effects  of  Exploding  Propellant  Tanks  and  Gas 
Storage  Vessels,"  NASA  CR-134906,  November  1975. 

1.4  X  10*  ft-lb  -  1  pound  TNT 

•  Kinney,  G.F.,  et  al,  "Explosive  Shocks  In  Air,"  1?85 
1.545  X  10*  ft-lb  "  1  pound  TNT 

•  U.S.  Air  Force,  "Air  Force  System  Command  Design  Handbook  1- 
6,"  December  1982. 

1.6  X  10*  ft-lb  -  1  pound  TNT 

As  shown  above,  there  are  many  conversion  factors.  T..e  average  of 
these  seven  is  1.477  x  10*  ft-lb  per  pound  TNT.  Kinney's  value  differs  by 
4.4  percent  from  the  average.  By  using  1.545  x  10*  ft-lb/pound  TNT 
(i.e.,  Kinney),  a  smaller  TNT  equivalence  is  obtained  and,  therefore, 
lower  overpressures.  Although  this  is  not  the  conservative  approacli,  it 
is  felt  that  being  able  to  relate  to  the  standard  definition  where  a  ton 
of  TNT  is  defined  as  the  energy  release  of  one  million  kilocalories  i-s  a 
good  approximation.  Further  vessel  bur.st  experimentation  should  .settle 
this  debate. 
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2.3.3  Comparison  of  Experimental  Data  with  TNT  Equivalence 

In  order  to  perform  a  comparison  between  Pittman's  experimental 
data  [10  &  11)  and  TNT  equivalence,  tcUbles  (Appendix  A)  were  developed 
using  Pittman's  data,  vessel  burst  energy  (isothermal  and  isentropic) , 
TMT  equivalence,  and  overpressures  from  Kingery  lata.  The  following  is 
the  sequence  for  performing  the  comparison; 

1)  The  Pittman  data  of  vessels,  vessel  volume,  vessel  burst  pressure, 
vessel  material,  and  burst  medium  were  input  on  the  tables. 

2)  The  vessel  burst  energy  from  the  expansion  equations  (isothermal 
and  isentropic)  was  calculated  and  input  on  the  tables. 

3)  The  TNT  equivalence  was  calculated  based  on  1.54  5  x  10‘  ft-lb  per 
pound  TNT  from  Kinney  [4]  and  input  to  the  table.  Pittman  had 
calculated  a  TNT  equivalence  for  his  vessels  but  used  an  earlier 
approximation  of  1.428  x  10‘  ft-lb  per  pound  TNT,  resulting  in  the 
differences  shown  in  the  tables. 

4)  The  TNT  equivalence  wes  converted  to  overpressure  data  at  distances 
equal  to  the  Pittman  data.  This  was  performed  using  the  table 
Kingery  (9)  developed  (see  table  A-I  in  Appendix  B)  .  Kingery' s 
table  uses  a  term  (X)  defined  as  the  scaled  distance  from  ground 
zero  in  units  of  ft/lb  TNT.  Taking  the  distances  from  Pittman 
and  the  cube  roots  of  the  TNT  eqjivalence,  (X.)  was  calculated. 
Entering  Kingery' s  table  A-I  with  (X),  a  corresponding  overpressure 
(Ap.)  in  psi  was  obtained.  These  values  of  Ap,  were  recorded  on  the 
tables . 

5)  From  this  data,  overpressure  versus  distance  curves  were  plotted 
and  compared  against  actual  Pittman  data.  Figures  2-9  and  2-lC 
show  representative  plots  with  the  remaining  plots  being  contained 
in  Appendix  C . 

Figures  2-9  and  2-10  plot  overpressure  versus  di-stance  for  real 
argon  (Figure  2-9)  and  ideal  argon  (Figure  2-10)  for  both  isentropic  and 
isothermal  expansion  of  the  gas.  The  same  9  data  points  from  Pittman 
(vessel  6)  appear  on  each  figure.  The  curves  as3”me  that  all  of  the 
energy  from  the  burst  goes  into  the  blast  wave.  Bo  .i  the  isothermal  and 
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isentropic  curves  overestimate  the  pressure  seen  at  a  given  distance 
from  the  burst.  The  isentropic  expansion  case  for  real  argon  predicts 
the  experimental  data  most  nearly  accurately. 

The  following  conclusions  can  be  drawn  from  Figures  2-9,  2-10,  and 
Appendix  C:  1)  Overpressure  vs.  distance  curve  for  a  pressure  vessel 

rupture  does  not  correspond  to  that  for  a  TNT  explosion  of  equivalent 
energy.  2)  Using  TNT  burst  data  from  hemispherical  charges  is  not 
accurate  for  use  with  cylindrical  vessels  because  a  blast  wave  from  a 
pressure  vessel  rupture  is  strongly  influenced  by  burst  geometry.  3) 
The  factor  controlling  the  magnitude  of  the  airblast  overpressure  frot. 
pressure  vessel  rupture  in  the  Pittman  tests  was  the  jetting  direction 
of  the  high  pressure  gases.  4)  No  allowance  was  made  in  the  calcula¬ 
tions  for  the  reduction  of  energy  available  to  the  blast  wave  resulting 
from  the  energy  expended  in  vessel  tearing  and  fragment  acceleration. 
Brown  [6]  states  that  anywhere  from  20  to  80  percent  of  the  burst  energy 
goes  into  the  blast  wave. 

Estimates  of  the  energy  available  to  the  blast  wave  based  on  theory 
and  experiment,  for  munitions,  are  reported  in  literature.  The  most 
common  estimate,  according  to  Bagelso  (12),  of  the  energy  avail6dble  to 
the  blast  wave  is  the  modified  Fano  formula  (Joint  Munitions 
Effectiveness  Manual,  1970) . 

0.6  +  0.4 

Weff  is  the  energy  available  for  blast.  C/M  i.s  the 
charge  weight -to-metal  weight  ratio.  The  charge 
weight  is  taken  as  the  equivalent  TNT  weight. 

These  different  estimates  for  the  amount  of  energy  available  to  the 
blast  wave  have  to  be  verified  with  additional  vessel  burst  data. 


Weff  ” 
W 


Where : 


t 
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2.4  Experimental  Approach  and  Computer  Codes 


2.4.1  Experimental  ^proach 

The  Department  of  the  Air  Force  la  concerned  with  the  integrity  of 
both  flight  and  ground  support  pressure  vessels  and  any  hazards  that 
result  from  their  use.  Additionally,  they  must  be  able  to  accurately 
predict  the  potential  damage  if  a  pressure  vessel  were  to  burst.  As 
shown  In  earlier  sections,  an  accurate  method  to  relate  vessel  burst 
energy  to  blast  wave  overpressure  does  not  exist .  Nor  is  there  enough 
existing  experimental  data  to  validate  computer  codes.  Therefore,  an 
experimental  approach  should  be  used  to  produce  empirical  data  on  which 
estimates  of  overpressure  dajiiage  can  be  based  or  validated. 

Pittman's  experimental  approach  flO  &  11)  for  bursting  pressure 
vessels  and  collecting  data  uses  a  mathodology  that  follow-up  studies 
performed  today,  could  use.  in  the  10  years  since  Pittman's 
experiments,  there  have  been  advances  in  technology  with  regards  to 
instrumentation  such  as  transducers  used  to  measure  blast  wave 
phenometia.  These  advances  in  technology  can  be  integrated  into 
Pittman's  methodology  to  upgrade  the  experimental  approach  to  today's 
standards.  However,  a  few  problems  still  exist  in  today's  technology 
with  regards  to  blast  wave  measurement.  These  include  calibration  of 
the  blast  gages  and  measurement  of  peak  overpressure.  Blast  gages  arc 
typically  piezoelectric-type  pressure  transducers  intended  for  dynamic 
situations.  These  may  be  capable  of  only  semi-static  use,  whereas 
static  transducer  calibrations  are  typically  the  roost  precise.  The 
extra  effort  for  dynamic  calibrations  may  be  deslreable  to  check  rise 
time  and  ringing  which  wo\»ld  affect  accuracy  in  measuring  peak  over¬ 
pressure.  The  most  precise  measurement  of  overpressure  reqfi'ires  ideally 
an  inertialess  pressure  sensing  element  with  zero  time  delay.  In 
contrast,  arrival  times  may  be  more  simply  and  precisely  measured. 
Compromises  can  be  made  between  the  desired  characteristics  and  what 
feasibly  can  be  obtained  without  degrading  the  experimental  results. 
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2.4.2  Computer  Codes 


Figures  2-11  and  2-12  are  representative  of  the  types  of  curves 
that  can  be  generated  by  computer  modeling  of  the  blast  wave  from  a 
bursting  pressure  vessel.  The  codes  used  to  genert'e  these  curves  are 
two  of  the  available  codes  that  have  been  used  to  model  blast  effects 
from  explooives  and  pressure  vessels.  The  two  codes  examined  by  this 
study  are  the  Wondy  V  Code  and  the  Sharp  Shock  Wundy  Code. 

The  Wondy  v  code  was  selected  because  Pittman  had  used  an  earlier 
version  of  the  code  in  his  first  study  [10] .  The  Sharp  Shock  Wundy  Code 
was  used  because  it  includes  a  modification  that  may  more  accurately 
predict  the  sharp  shock  peak  at  the  leading  edge  of  the  blast  wave. 
(This  is  evident  by  examining  the  uhape  of  the  two  curves  in  Figure  2- 
12,  and  comparing  ther'  to  theoretical  blast  wave  curves  or  to  actual 
test  curves  from  vessel  burst  tests.)  Differing  opinions  exist  as  to 
the  actual  shape  of  this  curve.  Some  contend  that  a  vessel  burst  blast 
wave  ij  characterized  by  o  sharp  rise  to  a  peak  presc  ^re  and  an 
exponential  decay  thereafter  see  Figure  2-la,  while  ethers  say  there  is 
a  less  rapid  rise  tc  a  less  distinct  peak  with  a  gradual  decay  to 
ambient,  see  eig'rre  2-lb.  This  also  will  have  to  be  resolved  by  more 
♦esf ng,  because  peak  pressure  is  a  major  concern  in  assessing  the 
eff  ects  of  bla.st  waves. 

Figure  2-11  shows  the  peak  pressure  that  would  be  observed  at  any 
radial  location  from  the  vessel  as  determined  using  the  Sharp  Shock 
Wundy  Code.  Another  way  to  explain  this  curve  would  be  to  connect  the 
peaks  in  Figure  2-2  and  a  curve  similar  to  Figure  2-11  would  be 
obtained.  Figure  2-12  .shows  peak  pressure  (overpressure)  versus 
distance  for  a  .specific  instance  in  time  after  the  vessel  burst,  as 
determined  using  both  ccmiputer  codes.  The  Wundy  Sharp  Shock  Code  was 
used  for  nitrogen  and  the  Wondy  V  Code  was  used  fo'  ,aiar-i.  These  curves 
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OVERPRESSURE  (psi) 


PRESSURE  (bar) 


1.5-1 


Figure  2-12.  Pressure  vs.  Scaled  Distance  for  Argon  and  Nitrogen  at  approximately 
1025  M  sec/cm.  (From  ATR,  1987) 
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war*  ^v«loped  for  a  vessel  which  is  similar  to  the  vessels  used  by 
Pittwan.  The  followiny  are  the  vessel  characteristics  used  in  the 
computer  code  example: 

Spherical 

Volume  -  6  ft* 

Fluid  -  gaseous  nitrogen  (GN2),  gaseous  argon  (GAr) 

Burst  pressure  -  8000  psig 

Fluid  temperature  -  17°C  (63*?) 

Figure  2-11  shows  a  very  high  overpressure  (>100  psi)  at  the  vessel 
surface  R/Ro  -  1,  which  rapidly  decays.  However,  it  is  evident  from  the 
graph  that  significant  overpressure  exists  farther  out  than  113  ft  (R/R„ 
■*  100) .  In  this  case,  "significant"  means  overpressure  that  could  break 
windows,  which  is  possible  with  fairly  low  overpressures.  If  a  minimum 
acceptaible  overpressure  of  0.5  psi  is  used,  then  the  separation  distance 
from  this  graph  would  be  about  170  ft  (R/R,  ■  150)  .  This  is 

considerably  less  than  the  minimum  criteria  of  1250  ft  required  by 

current  practice  [13] . 

An  additional  aspect  to  be  considered  with  these  curves  is  the 
assumptions  made  in  doing  the  computer  calculations.  First,  the  blast 
wave  was  ass'omed  to  be  one-dimensional  (1-D)  (i.e.,  the  energy  from  the 
blast  wave  dissipates  equally  in  a  circle  from  the  blast  center) . 
Second,  the  calculation  assumes  that  all  the  energy  available  in  the 
pressurized  gas  is  dissipated  through  the  blast  wave.  Neither  of  these 
assumptions  arc  completely  accurate  in  a  real  .situation.  These 

assumptions,  along  with  ways  of  accounting  for  them,  are  the  subjects  of 

the  following  paragraphs. 

T)»e  fir.st  assumption,  that  the  blast  is  1-D,  assumes  that  the 
vessel  breaks  up  into  two  equal  halves  with  the  blast  wave  being 
released  equally  in  a  circle;  see  Figure  2-13.  This  is  an 
underconservative  assumption,  because  most  vessels  fail  at  a  localized 
area  of  higli  stress  .  Because  vhe  initial  failure  would  occur  at  a 
specific  location,  a  large  portion  of  the  energy  stored  in  the  vessel 


^  '1..'  / 


\  / 


’’v  / 
/ 

/ 

/ 

# 

I 

I 

I 

t 

I 

I 

_  \ 

\ 

t 

\ 

\ 


^  . -  V 


/ 

\ 


One-Dirhensional  Blast 


Realistic  Situation 

Figure  2-13.  Onr  'dimensional  Blast  vs.  Realistic  Situation 
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would  be  released  through  the  initial  failure,  in  one  general  direction, 
and  not  evenly  in  all  directions;  see  Figure  2-13.  Pittir.an  made  the 
following  observation  based  on  his  1976  experiments  [8] : 

The  blast  field  from  pressure  vessel  rupture  is  strongly 
influenced  by  burst  geometry.  The  vessels  used  in  this 
investigation  burst  into  two  pieces  creating  strong  argon 
motion  along  a  line  defined  by  the  vessel's  center  and  the 
point  on  the  vessel  where  rupture  first  occurred.  Airblast 
overpressures  measured  along  the  line  of  this  jet  were  more 
than  a  factor  of  4  greater  than  those  measured  in  the 
opposite  direction. 

To  still  use  a  simple  1-D  Code,  it  would  be  necessary  to  muJtiply 
the  energy  available  for  the  1-D  Code  by  a  correction  factor  to 
compensate  for  the  directionalization  of  the  energy  release  in  a 
realistic  situation.  Different  types  of  vessels  may  require  different 
correction  factors.  To  determine  what  correction  factor  would  be 
necessary,  a  testing  prcgreim  would  have  to  be  conducted  which  would 
compare  computer  runs  to  burst  tests.  An  alternative  method  to  using 
such  a  correction  factor  would  be  to  use  a  more  complex  two-  or  three- 
dimensional  computer  code  to  predict  overpressure.  But  a  cwo-  or  three- 
dimensional  code  could  only  be  "run”  on  a  supercomputer.  This  option 
was  not  pursued. 

The  second  assumption,  that  all  the  energy  is  dissipated  in  the 
blast  wave,  is  an  overconservati ve  assiomption.  It  was  stated  earlier  in 
Section  2.1.2  and  again  in  2.3.3  that  only  20-80  percent  of  the  energy 
is  dissipated  through  the  blast  wave.  This  overconservai  ive  assum.ption 
may  or  may  not  compensate  for  the  underconservative  assumption  discussed 
above.  In  either  case,  predictions  must  be  compared  to  tests  of  various 
shapes,  sizes,  materials  and  pressures  of  vessels  before  an  accurate 
prediction  method  can  be  developed. 

2  5  Evaluation 


The  two  principal  models  used  to  predict  the  blast  wave  generated 
by  a  pressure  vessel  burst  employ  isentropic  expansion  and  isothermal 
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expansion,  respectively,  of  rho  contained  gas.  In  each  model,  the 
energy  rele.ased  by  the  vessel  burst  is  determined  by  calculating  the 
work  performed  by  the  expanding  gas.  All  of  this  work  is  assumed  to 
supply  energy  to  the  blast  wave.  The  expansion  of  the  contained  gas  is 
neither  isothermal  nor  isentropic,  but  the  isentropic  expansion 
approximation  appears  to  model  the  blast  v  ive  more  closely,  using 
Pittman's  data. 

Assuming  that  the  expanding  gas  is  ideal  predicts  greater 
overpressures  than  if  real  gas  equations  are  used.  rn  both  cases, 
however.  the  predicted  overpressures  exceed  those  obtained 
experimentally . 

It  is  convenient  to  equate  the  energy  released  by  a  vessel  burst  to 
that  produced,  by  an  explosion  of  a  given  amount  of  TNT.  This  is  a 
useful  comparison  since  there  is  a  much  larger  body  of  experimental  data 
for  TNT  explosions  than  for  pressure  vessel  bursts.  Unlike  TNT 
explosions,  however,  a  significant  amount  of  the  energy  produced  by  a 
vessel  burst  is  expended  in  breaking  the  vessel  and  in  accelerating 
fragments . 

There  is  a  limited  amount  of  available  data  to  compare  the  blast 
wave  methodologies  to  experimental  or  accidental  vessel  bursts. 
Additional  testing  is  required  to  validate  a  methodology. 

Two  one-dimensional  blast  wave  computer  codes  are  available  to 
model  the  overpressure  as  a  function  of  distance  from  the  point  of 
vessel  burst.  The  one-dimensiona J  codes  assume  that  the  blast  wave 
progresses  radially,  that  is,  the  overpressure  at  a  given  distance  is 
the  same  in  any  direction.  Since  vessels  can  be  expected  to  fail  at  a 
flaw,  a  jetting  effect  should  be  produced  with  much  higher  overpressure 
along  the  line  of  the  jet  than  in  other  directions. 
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SECTION  3 .  FRAGMEiiTATION 


The  purpose  of  this  section  is  to  review  the  current  methodologies 
availad>le  to  predict  the  properties  of  the  fragments  generated  by  the 
burst  of  a  high  pressure  pneumatic  vessel.  The  section  begins  witn  an 
introduction,  which  briefly  describes  some  of  the  most  important  studies 
on  pressure  vessel  fragmentation,  and  which  is  followed  by  a  discussion 
of  terminology  relevant  to  fragmentation. 

Most  analyses  to  predict  the  initial  velocities  of  fragments  are 
based  upon  a  work  by  Taylor  and  Price.  This  analysis  and  subsequent 
analyses,  which  generalize  Taylor-Price  or  adapt  it  to  special  cases, 
are  reviewed.  Each  analysis  is  evaluated  and  compared  to  experimental 
data,  where  applicable.  One  model  is  also  discussed  which  utilizes  a 
real  gas  instead  of  an  ideal  gas  which  the  other  models  use;  a 
comparison  between  real  and  ideal  gas  effects  on  fragment  velocity  is 
made . 


once  the  initial  velocity  of  a  fragment  is  determined,  it  is 
desirable  to  calculate  its  range.  This  is  accomplished  using  the 
established  ballistics  equations  of  physics;  computer  codes  are 
identified  which  determine  fragment  range,  given  the  initial  velocity. 
Some  upper  limits  on  fragment  range  are  identified  for  special  cases. 

Experimental  data  on  fragmentation  is  presented  and  evaluated.  The 
lack  of  available  data  on  the  distribution  of  fragments  (the  number 
produced  and  their  masses)  is  identified. 

A  discussion  of  a  probabilistic  determination  of  fragment 
generation  and  trajectories  is  addressed.  This  methodology  involves 
determining  probability  functions,  based  on  experimental  data,  for 
fragment  number,  size,  shape,  mass,  initia)  velocity  and  range.  A  lack 
of  available  experimental  data  currently  limits  the  use  of  this 
methodology . 
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3.1  Introduction 


The  failure  of  a  pressure  vessel  not  only  generates  a  blast  wave 
but  also  produces  fragments,  usually  with  very  high  velocities.  These 
fragments  constitute  a  significant  hazard  to  personnel,  systems, 
components  and  structures  in  the  vicinity.  A  study  of  the  hazards 
associated  with  the  failure  of  a  vessel  must  include  consideration  of 
the  range  and  potential  impact  energy  of  each  of  the  fragments  produced. 

Brown  [6]  and  Balter  [14]  have  both  written  summaries  of  the  current 
{circa  1984)  state  of  the  art  in  fragmentation. 

Determining  the  velocities  of  fragments  from  pressure  vessel 
ruptures  has  been  undertaken  by  several  authors.  Taylor  and  Price 
(1971)  developed  equations  for  velocities  of  the  fragments  for  a 
spherical  pressurized  gas  vessel  bursting  into  2  h; Ives .  Their  analysis 
was  expanded  by  ^essey  and  Kulesz  (1976)  to  allow  for  the  rupture  of 
either  a  sphere  or  a  cylinder  into  any  number  of  fragments.  Baker 
(1975)  developed  two  computer  codes,  SPHER  and  CYLIN,  based  upon  the 
Taylor-Price  analysis,  to  determine  the  fragment  velocities  for  a  sphere 
and  cylinder,  respectively,  bursting  into  a  given  number  of  fragments 
(which  must  be  assumed)  ;  the  computer  code  FRAG  2,  also  based  upon 
Taylor-Price,  determines  fragment  velocity  and  other  parameters  for  a 
vessel  bursting  into  two  equal  fragments.  Baker,  et  al,  (1978) 
developed  the  code  UNQL  to  calculate  fragment  velocities  when  a 
spherical  or  cylindrical  vessel  bursts  into  two  unequal  fragments, 

Jager  (1981,  derived  the  fragment  velocity  for  -uo  cases:  a  jet 
propelled  fragment  torn  off  of  a  moored  pressure  vessel  and  a  self 
propelled  fragment  such  as  a  rocketing  cylinder.  Baum  (1983)  modified 
the  work  of  Kulesz,  et  al,  (1979)  to  determine  upper  limits  on  fragment 
velocity  by  analyzing  gas  properties  and  the  effect  of  the  rarefaction 
wave  generated  within  the  gas  by  tne  motion  of  the  accelerating 
fragments . 

Bessey  (1974)  also  extended  the  Taylcr-Price  auolysis  to  spherical 
vessels  containing  liquid  propellants.  The  and  oxidizer  are 
assumed  to  come  into  contact  and  partially  del-  u  .  He  compares  liis 
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reaulta  to  those  of  Taylor-Price  for  pressurized  gases.  Similarly, 
Kulesz  (1978)  derived  equations  for  the  roclteting  effect  of  a  large 
portion  of  a  liquid  propellant  vessel,  and  developed  the  computer  code 
THRUST . 

Determination  of  the  range  of  a  fragment  from  its  initial  velocity 
is  well  understood.  Ballistics  equations  from  physics  can  be  used  to 
determine  fragment  trajectories;  these  require  )cncwledge  of  initial 
target  speed  and  elevation  angle  and  also  the  drag  (or  lift)  coefficient 
for  its  orientation  in  fJ.ight.  Two  computer  codes,  TRAJE  and  FRISB, 
were  developed  by  Baker,  et  al,  (1975  &  1978)  to  determine  the  ranges  of 
drag  type  and  lift  type  fragments,  respectively. 

Presently,  there  is  no  theory  to  predict  the  distribution  cr 
dispersion  of  fragments;  that  is,  the  number,  size  and  mass  of  fragments 
must  be  determined  empirically. 

3 . 2  Terminology  [27] 

One  of  the  primary  hazards  from  vessel  failure  is  due  to  the 
fragments  that  are  generated.  Fragments  can  be  accelerated  to  very  high 
speeds  and  threaten  nearby  personnel,  systems  and  components. 
Additionally,  fragments  can  penetrate  nearby  storage  vessels,  causing 
vessel (3)  to  fail  and  possibly  leading  to  detonation  of  liquid 
propellant . 

Primary  fragments  are  portions  of  the  vessel  or  its  attachments 
that  are  accelerated  due  to  the  internal  pressure  of  the  vessel  upon 
failure.  As  few  as  one  fragment  may  be  produced  in  the  failure  of  a 
ductile  vessel.  Vessels  which  fail  in  a  brittle  manner  or  which  fail 
due  to  the  partial  detonation  of  a  liquid  ^^ropellant  will  generally 
produce  more  fragments. 

Secondary  fragments,  or  appurtenance.s,  are  produced  due  to  the 
action  of  the  blast  wave  on  nearby  objects,  such  as  tools,  components 
and  parts  of  buildings.  An  important  secondary  fragment  is  bro)cen  glass 
which  is  a  hazard  to  personnel  and  can  also  penetrate  thin-skinned 
objects  nearby. 
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Fragments  will  follow  a  trajectory  until  they  impact  an  object  or 
strike  the  ground.  The  forces  which  act  upon  them  are  gravity  and  the 
fluid  dyn2unic  forces:  drag  and  lift.  The  drag  or  fluid  force  acting  on 
a  fragment  is  given  by: 

F  -  1/2  C  p  Av’ 

where  p  -  air  density 

V  -  fragment  speed 
A  -  surface  area 

and  C  is  either  the  drag  coefficient  Q  or  the  lift  coefficient  C^,  as 
appropriate.  The  lift  and  drag  coefficients  must  be  empiricaJly 
detexrmined  and  graphed  or  tabulated. 

Fragments  with  lift  coefficient  greater  than  drag  coefficient  are 
called  lift  type  fragments;  they  are  usually  disk  shaped.  Fragments 
with  drag  coefficient  greater  than  lift  coefficient  are  called  drag  type 
fragments.  They  are  usually  chunked  shaped,  that  is,  all  of  t.he 
dimensions  are  of  the  same  order  of  magnitude.  Trajectory  calculations 
are  performed  differently  for  lift  type  fragments  and  drag  type 
fragments . 

The  relative  hazard  that  a  fragment  presents  is  dependent  upon  both 
its  trajectory  parauneters  and  its  impact  effects.  A  fragment  can  cause 
damage  by  either  penetrating  or  rebounding  from  an  object.  If  the 
fragment  passes  entirely  through  the  object,  the  process  is  called 
perforation.  If  the  fragment  fails  to  pass  entirely  through  the  object, 
but  displaces  a  portion  of  it,  it  is  termed  penetration.  Spalling  is 
the  process  by  which  the  impact  induces  compression  waves  in  the  object, 
causing  a  tension  failure. 

3 . 3  Determination  of  Initial  Fragment  Velocity 
3.3.1  Taylor  -  Price  Ana..j^.s 

The  velocities  of  2  fragments  generated  by  the  rupture  of  a 
spherical  pressure  vessel  were  predicted  by  Taylor  and  Price  [15]  in  a 
1971  .kSME  paper.  Their  work  was  based  upon  that  of  Grodzovskii  and 
Kukanov  (1965)  but  removed  two  re.st  rict  ions :  (1)  the  speed  of  the  gas, 
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•nd  therefore,  the  fragments  are  small,  relative  to  the  sonic  escape 
velocity  and  (2)  the  volume  between  the  separating  fragments  as  they  are 
accelerating  does  not  change  significantly  from  the  original  volume. 

The  Taylor-Price  analysis  assumes  that  the  two  fragments  move  in 
opposite  directions  in  a  vacuum  and  that  an  ideal  gas  escapes 
perpendicularly  to  the  motion  of  the  fragments.  The  velocities  are 
determined  for  the  fragments  for  two  cases:  isothermal  and  adiabatic. 

The  isothermal  c-se  yields  higher  velocities  for  the  fragments  than 
does  the  adiabatic  case.  This  is  not  surprising  since  the  isothermal 
case  predicts  higher  available  energy  from  the  blast .  Both  cases 
predict  lower  velocities  than  those  obtained  by  Grodzovskii  and  Kukanov. 

In  the  isothermal  case,  they  determined  that  as  the  dimensionless 
mass  a  ”  P„  Vq  =•  C„  becomes  large,  the  maximum  possible  fragment 
M.a,’  KM, 

velocity  v^  approaches  a  constant :  v^  -  a„ 

where  a,  “  sonic  velocity  of  gas 
C„  »  mass  of  gas 
K  “  ratio  of  specific  heats 
M,  -  total  mass  of  vessel  and  gas 

This  maximum  fragment  velocity  corresponds  to  the  maximum  velocity 
attainable  in  nozzle  expansion.  Note  that  Ct  is  large  in  highly 
stressed,  lightweight  vessels. 

The  analysis  by  Taylor  and  Price  has  been  used,  extended  or 
modified  by  several  other  authors  in  order  to  inprove  on  the  assumptions 
or  to  adapt  the  analysis  to  other  cases. 

One  such  extension  of  Taylor-Price  was  performed  by  Baker,  ot  al, 
(1975),  and  covers  cylinders  as  well  as  spheres. 
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3.3.2  Modifications  to  Taylor  -  Price 


Baker  [5]  assumed  that,  in  the  case  of  a  sphere,  it  fra9ments  into 
n  fragments  of  circular  projection  which  travel  radially  without 
tumbling.  A  cylinder  shell  is  assumed  to  fragment  into  n  lengthwise 
strips  which  travel  radially;  motion  of  the  two  heads  is  not  considered. 
The  energy  of  the  contained  gas  is  partitioned  between  the  kinetic 
energy  of  the  fragments,  the  energy  of  the  escaping  gas  and  the  energy 
of  expansion  of  the  gas.  Strain  energy  in  the  vessel  walls  prior  to 
burst  is  neglected;  however,  for  steel  vessels  containing  gas  at 
pressures  of  several  thousand  psi,  this  energy  is  on  the  order  of  0.1%  - 
0.2%  of  the  gas  energy. 

Baker  developed  the  computer  codes  SPHER  and  CYLIU  to  solve  for  the 
nondiroensional  velocities  of  the  fragments  for  the  spherical  and 
cylindrical  vessels,  respectively.  Both  codes  use  the  Runge-Kutta 
integration  technique  to  solve  simultaneous  nonlinear  differential 
equations  for  the  velocities.  The  number  of  fragments  generated  by  the 
burst  must  be  assumed  and  used  as  an  input  parameter.  Baker  used  these 
Codes  to  solve  for  the  velocities  in  a  larce  number  of  cases.  These 
solutions  were  used  to  generate  curves  which  could  be  employed  to 
determine  fragment  velocities. 

The  solutions  demonstrated  that  as  the  number  of  fragments  produced 
increases,  the  maximum  fragment  velocity  increases  to  a  maximum;  this 
occurs  when  10-30  fragments  are  produced. 

The  SPHER  code  was  run,  and  the  curves  were  also  used,  to  calculate 
fragment  velocities  tor  4  cases  which  were  measured  by  Pittman  [10] 
(discussed  in  Section  S.*!)  .  The  SPHER  code  and  curves  underestimate  the 
velocity  of  the  smaller  sphere  by  about  10%.  They  underestimate  the 
velocity  of  the  larger  sphere  by  20-25%  as  measured  by  the  breakv/ire 
system  but  are  very  close  as  measured  by  the  strobe  photographic  system. 
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The  measured  and  calculated  velocities  are: 


Sphere  Characteriatica  Pressurizing  Gas  Initial  Fragment  Velocities 
Wall 


Radius 

Thickness 

Pressure 

V (Pittman) 

V(Code) 

V  (Carves) 

cm 

cm 

Tjjge 

Pa 

m/s 

m/s 

m/s 

11.7 

0.274 

V, 

5.51x10' 

366+15 

352 

338 

34.3 

0.919 

«, 

5.51x10’ 

342+30' 

339 

322 

34.3 

0.919 

N. 

5.51x10’ 

426+27’ 

339 

322 

34.3 

0.919 

N, 

5 . 51x10’ 

448+30' 

339 

322 

‘This  value  was  based  on  velocity 
photographic  technique. 

measurements 

using 

a  strobe 

‘These  values  were  based 

measurement  techniques. 

on  velocity 

measurements 

using 

brea)twire 

Baker  also  modified  the  Taylor-Price  method  to  solve  for  fragment 
velocities  when  a  vessel,  sphere  or  cylinder,  bursts  into  two  equal 
halves.  The  following  assumptions  were  made: 

(1)  Vessel  breaks  into  two  equal  halves  along  a  plane 
perpendicular  to  cylindrical  axis  and  the  fragments  are 
driven  in  opposite  directions. 

(2)  Contained  gas  obeys  the  Ideal  gas  laws. 

(3)  Escaping  gas  travels  perpendicular  to  the  direction  of  motion 
of  the  fragments  with  local  sonic  velocity,  and  enters  a 
vacuum. 

(4)  Energy  necessary  to  break  rhe  vessel  walls  is  negligible 
compared  to  the  total  system  energy. 

(5)  Drag  and  lift  fortes  are  negligible  during  the  time  period  in 
which  the  fragment.^  attain  their  maximum  velocities. 

A  computer  code,  FRAG  2,  was  developed  to  solve  for  veloc  ty, 
acceleration,  distance,  time  and  ambient  pressure.  The  program  was  run, 
varying  input  para-meters,  for  the  following  conditions: 

(1)  Vessal  is  made  of  titanium  or  titanium  alloy. 

(2)  Vessel  walls  are  of  uniform  thickne.as. 

(3:  Vessel  has  hemispherical  heads. 
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The  following  conclusions  can  be  drawn: 


(1)  As  pressure  increases,  fragment  velocity  increases  at  a 
decreasing  rate. 

(2)  As  gas  density  increases,  sonic  velocity  increases,  fragment 
velocity  increases,  and  the  rate  of  change  of  velocity  with 
respect  to  pressure  (dv/dP)  decreases. 

(3)  As  vessel  thickness  to  diameter  decreases  (i.e.,  for  thinner 
wall  vessels)  ,  fragme-tit  velocity  increases  and  the  rate  of 
change  of  velocity  with  respect  to  pressure  (dv/dP) 
decreases . 

(4)  As  vessel  length  to  thickness  decreases,  fragment  velocity 
increases,  particularly  at  lower  pre.ssures. 

Gases  that  were  chosen  for  the  calculations  were  air,  hydrogen, 
xenon  and  carbon  dioxide. 


3.3.3  Self-Propelled  and  Jet-Propelled  Fragments 

Jager  [16,  17]  derived  corre  relatively  simple  equations  to 
determine  the  maximu,::  velocity  that  a  fragment  migliv  obtain  for  two 
cases.  The  first  case  covers  self  propelled  fragments  such  as  a 
rocketing  gas  cylinder;  the  second  case  addresses  jet  propelled 
fragir.ents  such  as  a  flange  torn  off  of  a  moored  vessel. 


The  analysis  ignores  gravity  over  the  time  period  during  which  the 
fragment  is  accelerated  to  its  final  velocity  It  is  assumed  that  both 
the  choked  flow  velocity  and  density  --.ecrease  exponentially  as  the 
vessel  discharges  with  time. 


The  maximum  velocity  attainable  by  a  self  propelled  fragment 


giver,  by 
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where  m,  -  mass  of  gas  initially  in  vessel 
m,  -  mass  of  rocketing  f.agment 
K  -  ratio  of  specific  heats 
a,  -  gas  sonic  speed  (initial) 

P„  “  ambient  pressure 

Pi  -  initial  pressure  in  vessel 

3.3.4  Limiting  Case  Energy  Considerations 

Baum  [18]  (1984)  considers  limiting  cases  to  determine  ragment 

velocities.  Tf  a  gas  at  initia'  pressure  P„  and  initial  volume  v„ 
expands  adiadDaticaliy ,  the  maximum  expansion  work  E  is  given  by 

E  -  P<,_V„ 

K-1 

where  K  is  the  ratio  of  the  specific  heats. 

In  the  case  of  an  ideal  gas  (i.e.,  initial  pressure  is  sufficiently 
low  that  molecular  forces  and  finite  molecular  size  do  not  influence 
behavior  of  the  gas,  and  condensation  does  not  occur  at  the  end  of 
expansion)  ,  the  maximum  proportion  k  of  thi.s  energy  which  could  appear 
as  kinetic  energy  is 

where  P,  is  the  pressure  external  to  the  vessel . 

The  me  rmum  kinetic  energy  available  for  fragments  is  obtained  by 
multiplying  ttie  proportion  k  by  the  available  energy: 

-  P._V.  k 
K-1 

Baum  [21]  (1988)  has  also  performed  burst  tests  for  a  number  of 
different  cases  oi  fragment  types  and  vessel  contents  for  both  spherical 
and  cylindrical  vessels.  Initial  fragment  velocities  were  measured  and 
evaluated.  In  general,  initial  fragment  velocity  is  a  function  uf : 


•  velocity  of  sound  in  the  gas 

•  initial  fragrrient  acceleration 

•  ratio  of  initial  internal  vessel  pressure  and  external 

pressure 

•  vessel  lenqt.h  to  radius  ratio 

•  isentropic  expansion  coefficient  of  gas 

Based  upon  results  from  his  own  experiments  and  those  performed  by 
others,  Bavuti  recommends  upper  limits  on  fragment  velocitie  j  for  se'^eral 
cases.  Cases  for  vessels  containing  ideal  gases  are  summarized  beljw; 


Vessel 

Geometry 

cylindrical 

cylindrical 

cylindrical 


Failure 

Type* 

D  or  D/B 

D  or  D/B 

D/B  or  B 


spherical  B 


cylindrical  D 


Fragment 

Type 


Upper  L.imit  Velocity*  * 


end  cap 


rocket 


multiple  fragments 
from  vessel  dis¬ 
integration 


2.18 


1/2 

2  a.(F) 


0.55 

0.88  a,  (F) 


0.55 

multiple  fragments  0.88  (F) 

from  vessel  dis- 
integi ation 

1/2 

whole  vessel  with  0.17  (2E/m) 

axial  split 


Notes : 

(1)  Failure  type  is  ductile  (D) ,  brittle  (B)  or  ductile/brittle  (D/B) 

(2)  F  is  the  dimensionless  initial  fragment  acceleration  given  by 
1  -  P,  A  R 

ma,*  where 

P„  •*  initial  pressure 
A  -  fragment  projected  area 
l;  -  vessel  radius 
m  -  fragment  mass 
a,  -  velocity  of  sound  in  gas 
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(3)  other  teen's 

L  -  vessel  length 
E  ••  maximum  expansion  work 


3.3.5  Real  Gas  Effects 

Hledermann  [19]  employed  a  quasi-steady  expanding  cavity  model  similar 
to  Baum  to  evaluate  the  effects  that  a  real  gas  has  upon  fragments.  Each 
of  the  preceding  models  assximed  that  the  gas  in  the  vessel  was  ideal.  This 
is  a  good  assumption  for  low  pressure  gases.  At  higher  pressures,  real  gas 
effects  are  more  significant.  Where  Baum  used  the  ideal  gas  law, 
Wiedermann  employs  the  Nobel-Abel  and  Van  der  Waals  equations  of  state. 

Wiedermann  developed  non-dimensional  equations  which  he  solved  for 
nor.-dimensional  fragment  velocity.  The  equations  were  solved  for  a  range 
of  values  for  the  following  dimensionless  parameters: 

acceleration, 

co-volume  (difference  in  volume  between  real  gas  and  its  associate^ 

ideal  gas) , 

ratio  of  specific  heats  (adiabatic  exponent),  and 

gas  volume  factor. 

The  results  indicated  a  significant  reduction  in  fragment  velocity 
for  real  gas  in  the  bursting  vessel  as  opposed  to  an  ideal  gas. 

3 . 4  Determination  of  Fragment  Range 
3.4.1  Theory 

Once  the  initial  velocity  of  a  fragment  h  been  determined,  its 
range  may  be  determ-’ ned  through  ballistics  calculations.  The  equations 
of  motion  for  a  fragment  can  be  written  in  x  and  y  coordinates  where  x  is 
the  horizontal  range  and  y  is  the  altitude. 


5T 


For  a  drag  type  fragment,  the  acceleration  that  it  experiences  is 
due  to  gravity  (in  the  -y  direction)  and  drag  (opposing  the  velocity) . 


F,  -  X  -  -  Cp^_p  v'  cos  6 

m  2m 


sin  6  -  g 


where  m  mass  of  fragment 


X  -  horizontal  acceleration 


y  -  vertical  acceleration 
Cj,  -  drag  coefficient 

Aq  -  area  exposed  to  drag  or  reference  area 

p  -  air  density 

v  ■  velocity 

0  -  trajectory  angle 

g  -  acceleration  due  to  gravity 

F  -■  net  force  acting  on  the  fragment 


In  the  case  of  a  lift  type  fragment,  Co  Ao  are  replaced  by 

and  Ao,  the  lift  coefficient  and  area  exposed  to  lirt,  respectively.  If 
lift  and  drag  are  of  the  same  order  of  magnitude,  both  terms  might  be 
used  in  the  equations. 


The  initial  conditions  (boundary  values)  are  as  follows.  At  time 
t  -  0  :  X  (0)  -  V,  cos  0  „ 

y(0)-  V,  sin  0  „ 

where  x  -  horizontal  velocity 

Y  -  vertical  velocity 

v„  -  initial  velocity 

0„  -  initial  trajectory  angle 


The  velocity  v  may  be  related  tc  the  horizontal  velocity  x  and 
the  vertical  velocity  y  by 

I 

v  -  (x’  'f) 

The  trajectory  of  the  fragment  is  determined  by  solving  the  two 
differentia.  c;piatic:;3  cimu]  taneouoly .  This  is  generally  done  thrc';gh 
the  use  of  a  computer  code. 
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The  solutions  will  generally  give  velocity  and  horizontal  and 

vertical  distances  as  functions  of  time.  The  range  would  be  the 

horizontal  distance  when  the  vertical  distance  is  0. 

3.4.2  Computer  Codes  for  Fragment  Range 

Kulesz,  et  al  [5]  developed  the  computer  codes  FRISB  and  TRAJE  to 
determine  the  ranges  of  fragments.  FRISB  is  used  for  disc  shaped 

fragments  with  diameter  at  least  5  times  greater  than  thickness;  lift 
effects  are  considered  by  FRISB.  TRAJE  is  used  for  chunky  sha’-ed 
fragments  with  all  three  dimensions  of  the  same  order  of  magnitude. 

Lift  effects  are  ignored  in  TRAJE.  The  FRISB  code  can  also  be  used  for 

chunky  shaped  objects  by  setting  the  lift  coefficient  equal  to  zero. 

FRISB  and  TRAJE  both  employ  the  Runge-Kutta  method  to 
simultaneously  solve  two  second  order  differential  equations  for 
fragment  velocities.  The  velocities  are  numerically  integrated  in  order 
to  determine  the  range. 

The  TRAJE  code  can  access  a  subroutine  that  allows  for  rotor 

effects  such  as  a  helicopter  blade  would  demonstrate. 

Kulesz  ran  the  code  TRAJE  for  different  values  of  the  input 

variables : 

initial  trajectory  angle, 
initial  velocity, 
mass, 
area,  and 

aspect  ratio  (diameter/thickness) 

The  computer  run  results  were  used  to  generate  curves  that  allow  for 
graphical  determination  of  range.  Sample  checks  were  run  by 

interpolating  results  from  curves  and  comparing  these  results  to  the 

TRAJE  predictions.  The  error  ranged  from  5%  to  80%.  The  greatest 
errors  occurred  for  low  ratios  of  fragment  mass/area. 
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3.4.3  Upper  Limits  on  Range 


Baum  (18]  (1984)  recommends  upper  limits  on  range  similar  to  the 
upper  limits  for  initial  fragment  velocity.  In  the  case  where  the  drag 
force  is  less  than  the  fragment's  weight  and  t  e  lift  force  is 
negligible,  an  upper  limit  on  range  can  be  determined  by  solving  the 
ballistics  equations  by  setting  drag  and  lift  coefficients  equal  to  zero 
and  assuming  the  optimal  initial  trajectory  angle  of  45°.  Then  the  range 
is  given  by: 


g 


For  non-tumbling  fragments  with  significant  lift  force  but  drag 
force  still  less  than  weight,  the  upper  limit  on  range  is  given  by: 

X  -  3__^' 

g 


The  upper  limit  on  range  for  a  large  fragment  from  the  ductile 
failure  of  a  vessel  can  be  determined  by  assuming  that  the  )cinetic 
energy  of  the  fragment  is  20%  of  the  maximum  expansion  wor)c  available 
(E)  and  that  the  fragment's  mass  is  one  tenth  of  the  mass  (M)  of  the 
vessel . 


X  -  4  E 
Mg 

The  upper  limit  on  range  for  any  of  the  small  fragments  generated 
by  the  brittle  failure  of  a  vessel  can  be  determined  by  assuming  that 
the  total  )tinetic  energy  of  the  fragments  is  40%  of  the  maximum 
expansion  wor)c  and  that  each  of  the  fragments  has  the  same  initial 
velocity . 


X  -  0 .8  E 
Mg 
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3 . 5  Experimental  Data 
3.5.1  Pittman,  1972 

In  1972,  Pittman  [10]  burst  five  vessels  of  titani'um  alloy  (Ti-6A1- 
4V)  ;  three  of  these  were  spheres  used  to  store  high  pressure  gas  while 
two  were  cylinders  used  to  store  liquid  propellant.  Design  burst 
pressures  were  7500  psig  or  8000  psig  for  the  spherical  vessels  and  460 
psig  for  the  cylindrical  vessels.  All  five  vessels  were  pressurized 
with  gaseous  nitrogen  to  burst . 

The  fragment  recovery  system  utilized  thirty  6-inch  thick  wallboard 
panels  to  stop  the  fragments.  Two  fragment  velocity  measurement  systems 
were  used.  The  first  system  used  breakwires  in  the  panels  to  trip 
electric  counters.  The  second  system  used  two  cameras  with  strobes. 
Speed  was  estimated  by  computing  the  difference  in  position  for  each 
identifiable  fragment  between  two  stroboscopic  flashes  assuming  straight 
line  fragment  travel.  At  least  one  of  the  velocity  measurement  systems 
functioned  for  each  vessel  burst. 

The  number  of  fragments  that  each  vessel  burst  into  was  estimated 
by  calculating  the  percentage  of  the  vessel  weight  that  was  recovered 
and  taking  the  ratio  of  the  number  of  fragments  recovered  to  the  total 
vessel  weight.  Results  for  the  vessels  are  as  follows: 


Vessel 

Usual 

Content 

Number  o: 

A 

Ti.  f"  ^ 

A- 

61 

B 

Oxidizer 

57 

C 

GHE 

38 

D 

GHE 

47 

E 

GHE 

48 

61 


upon  detailed  analysis  of  the  fragments  produced,  Pittman 
concluded: 


•  The  method  of  estimating  the  number  of  fragments  is 
incorrect . 

•  The  fragment  distribution  is  not  homogeneous  over  a  solid 
angle  described  by  a  spherical  shell. 

•  Fragments  found  within  the  arena  went  straight  up  and  fell 
back . 

Measured  fragment  velocities  are  recorded  in  Table  3-1 
3.5.2  Pittman,  1976 

Pittman  (11]  in  1976  pressurized  seven  spherical  vessels 
constructed  of  T-1  steel  with  argon  until  each  failed.  Each  had  a 
capacity  of  1  cubic  foot.  The  design  burst  pressures  were  15  /.si,  30 
ksi  and  50  ksi.  Each  vessel  burst  into  2  fragments.  Of  the  14 
fragments,  10  were  recovered.  Data  for  the  remaining  4  fragments  were 
estimated.  The  weights  of  the  fragments  varied  from  51  pounds  (from  a 
15  ksi  vessel)  to  271  pounds  (from  a  50  ksi  vessel) . 

The  fragment  recovery  and  velocity  measurement  system  was  designed 
assuming  that  the  vessels  would  fail  in  a  brittle  manner,  and  produce 
several  fragnients.  Velocity  screens  were  set  up  to  stop  the  fragments. 
Breakwires  in  the  screen  would  trip  electric  counters  upon  impact, 
allowing  velocity  to  be  determined  from  the  arrival  time. 
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Strobe  S 

ystem  Velocities 

Tan)c 

Brer.kwire 

Fragment 

Time 

Distance 

Velocity 

Velocity 

Number 

Interval  (m3) 

Traveled  (ft) 

VA  * 

0  to  1.50 

1.52 

1010 

1.50  to  3.00 

1.56 

1040 

1 

3.00  to  5.00 

1.50 

750 

0  to  5.00 

4 . 58 

918 

A 

No 

0  to  1.50 

1 . 56 

1040 

Data 

1.50  to  3.00 

1.61 

1070 

2 

3.00  to  5.00 

1.86 

930 

0  to  5.00 

5.03 

1000 

0  to  3.00 

3.10 

1030 

3 

3.00  to  5.00 

1.73 

865 

0  to  5.00 

4.83 

967 

0  to  3.00 

2 .76 

920 

1 

3.00  to  4.50 

1 .29 

860 

B 

1215 

0  to  4.50 

4.05 

900 

+  50 

0  to  3.00 

2.96 

986 

2 

3.00  to  4.50 

1.32 

880 

0  to  4 . 50 

4.28 

950 

C 

1270 

1 

0  to  3.00 

3.60 

1200 

+80 

3.00  to  4.50 

1.53 

1020 

0  to  4.50 

5.20 

1150 

D 

1400 

No  Strobe  Data 

+  90 

E 

1470 

+  100 

No  Strobe  Data 

*  VA  -  Averagt!  velocity  in  ft/sec  over  the  time  interval  given 


Table  3-1.  Pittman  1972  Fragment  Data  [10] 
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ro  N)  N> 


All  of  the  veasels  failed  in  a  ductile  manner  and  burst  into  2 
fragments.  Only  1  of  the  14  fragments  tripped  the  breakwire  to  record  a 
velocity.  Velocities  of  the  other  fragments  were  measured  based  on  a 
pressure-temperature  history. 


SHOT 

# 


1 

1 


3 

3 

3 

3 

4 
4 
4 
4 


6 

6 

7 

7 

7 


Fragment  weights,  measured  velocities  and  ranges  are  as  follows. 


VESSEL  DESIGN 
BURST  PRESSURE 
(PSD 

FRAGMENT 

WEIGHT 

(POUNDS) 

FRAGMENT 

VELOCITY 

(FEET/SEC) 

DISTANCE 

FRAGMENT 

TRAVELED 

(FEET) 

15, 000 

52 . 4 

310'*' 

52.4 

321'*’  BOTTOM 

HALF 

665‘” 

TO?  HALF  NOT  FOUND 

1680"' 

15, 000 

52 .75 

324"’ 

52.75 

330'*’  BOTTOM 

HALF 

317'*’ 

552'*’ 

TOP  HALF  NOT  FOUND 

1380"’ 

30, 000 

147.5 

354"’ 

147.5 

373"’  BOTTOM 

HALF 

INSIDE  THE 

147.5 

348'*’ 

ARENA 

147 . 5 

338"’ 

TOP  HALF  NOT  FOUND 

1136"’ 

50, 000 

258 

215'*’ 

258 

216"’ 

INSIDE  THE 

258 

210'*’  LARGE  PORTION 

ARENA 

258 

221"’ 

SMALL  PORTION  NOT  FOUND 

722"’ 

50, 000 

271 

270'*’ 

INSIDE  THE 

142 

ARENA 

15, 000 

51.1 

NONE 

51 . 1 

NONE 

30,000 

108 

290'*’ 

INSIDE  THE 

108 

250'*’ 

ARENA 

93.5 

NONE 

Notes : 

(1)  VELOCITY  MEASUREMENTS  MADE  FROM  P-T  HISTORIES 

(2)  VELOCITY  MEASUREMENTS  MADE  FROM  BREAKWIRE  SYSTEM 

(3)  MEASURED  DISTANCE  TRAVELED 

(4)  ESTIMATED  DISTANCE  TRAVELED 


Table  3-2.  Pittman  1976  Fragment  Data 
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3.5.3  Jager,  1981 


Jager  (17,20)  in  1981  measured  the  velocities  of  single  fragments  jet 
propelled  from  26  .steel  compres  ed  air  bottles.  Prior  to  the  experiments, 
contours  of  the  desired  fragments  were  machined  into  each  bottle.  The 
fragments  were  either  5cm  by  5cm  squares  or  15cm  by  15cm  squares.  Masses 
ranged  from  0.16)cg  to  150kg.  In  some  cases,  however,  the  vessel  was  cut 
in  half  in  order  to  obtain  a  self  propelled  fragment.  The  bottles  were 
r nchored  in  place  to  minimize  movement . 

Velocities  were  measured  by  photographing  the  fragments  as  they  moved 
past  scaled  distance  marks.  The  fragments  moved  past  the  scale  in  15  of 
26  tests. 

Fragment  velocities  ranged  from  20  m/s  to  280  m/s.  The  table 
following  lists  both  measured  velocities  and  velocities  calculated  using 
equations  developed  by  Jager  (discussed  in  Section  3.3.3). 

3 . 6  Fragment  Distribution 

There  are  currently  no  theories  to  predict  the  number,  sizes  and 
masses  of  fragments  produced  as  a  result  of  a  pressure  vessel  burst . 
Studies  in  this  area  have  centered  on  statistically  curve  fitting  data 
from  the  experiments  by  Pittman. 

Baker,  et  al  (5)  analyzed  data  Pittman  [10]  obtained  by  bursting 
titanium  vessels.  They  found  that  the  logarithms  of  the  fragment  masses 
follow  a  normal,  or  Gaussian,  distribution.  They  plotted  two  mass 

distribution  curves - one  for  the  percentage  of  fragments  with  a  weight  < 

W  versus  weight  W  and  the  other  for  the  mean  fragment  weight  versus 
normalized  yield.  Normalized  yield  -  PV/E„  where  P  and  V  are  vessel 
pressure  and  volume,  respectively,  and  E„  is  the  energy  of  detonation  of  1 
greim  of  TNT  (4190  joules)  .  There  were  insufficient  data  points  on  the 
latter  curve  to  determine  if  it  was  a  straight  line  on  a  log-log  scale. 
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Tes 

no 


(bars) 


F.ICO 

(mM 


RE^iARKS 


_L- 


Vp  (m/s) 
calc,  imeas. 


JET  PROPELLED  FRAGMENTS 


JcrJ; 


Crr‘ 


1 

2 

3 

H 

116 

1 1 7 
215 

50 

50 

50 

5.8 

5.9 
10.7 

1 .43 
1.38 
1.35 

2.5 

2.5 

2.5 

205 

206 
226 

vessel  burst  in  3  pcs. 

rrr-- 

B 

KM 

116 

50 

5.8 

0.16 

0.35 

178 

impact  before  meas. 

174 

50 

S.7 

0.16 

0.35 

267 

158 

not  clearly  identified 

B 

mM 

381 

50 

19.0 

0.15 

0.35 

277 

280 

burst  into  2  pcs. 

- 

)  ’•■■iszs'. 

hot  torr 

(pure  lu 

jet-prspelied,  no  pluo--: 

upe  phsse) 

7 

100 

■  i| 

SO 

5.9 

12.0 

3.5 

97 

84 

no  rigid  support 

Bil 

SO 

11.7 

12.66 

3.5 

129 

119 

no  rig)(J  support 

Wm 

28  5. 

so 

16.5 

12.49 

3.5 

148 

125 

no  rigid  support 

B 

25 

30 

220 

6.6 

21 .8 

18.4 

96 

90 

rigid  support 

rocketing 

FRAGMENTS 

1 

100 

127 

50 

3.1 

33.2 

3.5 

57 

53 

1/2  vessel 

12 

29d 

340 

50 

8.2 

31. C 

3.5 

161 

156 

1/2  vessel 

w 

100 

119 

50 

2.0 

25.1 

3.5 

50 

60 

1/.3  vessel 

1  ^ 

30 

6.2 

119.5 

13.4 

30 

26 

whole  vessel  without 

bOttOTl 

15 

25 

31 

220 

6.8 

1 56 

>30 

24 

16 

vessel  cut  open  length- 

wise,  flying  side- 

_ 1 

wards  and  rotating 

•Flying  third  of  vessel  treated  as 

jet-propelled  from  remaining  2/3-vessel 

(with  (11-  =  3 

.9  Kg)  would  give  calculated  v 

59  m/s. 

1  Jw-ro 

P,  - 

pressure  of  gas 

in  vessel .  o,  - 

density  of  gas 

in  ves 

sel  ,  V,  -  volume 

1  capacity  of  vessel 

,  m  -  mass  of  gas  in  vessel  or 

vessel 

portion,  m  -  mass  of 

1  fragment, 

1 _ _ _ 

Vp  -  peak  velocity  of  fragment. 

J 

fl7] 

Table  3-3.  Jager  Fragment  Data 
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3 . 7  Probabiliatic  Methodology 

Rather  tnan  base  a  study  of  fragmentation  cn  analytical  methods  and  a 
"worst  case"  scenario,  a  probaibil^ stic  methodology  might  be  employed.  The 
probabilistic  methodology  relies  upon  empirical  data  in  order  to  generate 
probabilistic  den.sities  or  curves  for  damage  sustained  as  a  consecpience  of 
a  bursting  vessel's  fragments.  The  validity  of  a  probabilistic  assessment 
is  a  function  of  the  availability  and  accuracy  of  the  data  and  also  the 
choice  for  an  acceptable  level  of  damage  to  be  sustained. 

A  prob2d3ilistic  assessment  for  damage  as  a  result  of  vessel 
fragmentation  has  been  performed  by  Sundararajan  and  Roolcer  (23)  (1984)  . 
Their  approach  divides  the  assessment  into  two  phases: 

1.  generation  of  fragments 

2.  fragment  trajectories 

In  the  first  phase,  the  failure  probedjility  of  a  pressure  vessel  is 
examined  for  a  hypothetical  failure,  the  distribution  of  the  number  of 
fragments  and  their  weights,  sizes  and  shapes,  and  fragment  velocities 
(including  ejection  angles)  are  analyzed  to  develop  probability  functions. 

The  second  phase,  determination  of  fragment  trajectories,  could  be 
accomplished  through  one  of  two  approaches:  a  Monte  Carlo  simulation 
method  or  a  semi-analytical  method.  The  Monte  Carlo  simulation  method 
consists  of  performing  a  series  of  trials  to  determine  the  ranges  of 
fragments  by  sampling  input  fragment  speeds  and  ejection  angles  according 
to  the  probability  density  functions  obtained  in  the  first  phase.  From 
the  series  of  trials,  a  strike  probability  contour  map  is  drawn;  this 
contour  map  allows  determinations  to  be  made  for  the  relative  hazards  of 
locating  personnel  or  equipment  around  the  vessel. 

The  semi-analytical  method,  alternatively,  might  be  used  to  determine 
probability  density  functions  of  impact  velocity,  fragment  orientation  and 
weight  for  targets  located  near  the  vessel. 
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A  probabilistic  assessment  has  been  performed  by  Huang  ’,23]  to 
determine  the  impact  effect  of  debris  from  a  building  due  to  the 
detonation  of  stored  munitions.  Since  much  data  is  available  on  munitions 
blasts,  it  is  possible  to  develop  a  probabilistic  assessment  whereas  for 
vessel  bursts,  insufficient  data  is  available.  This  methodology  might  be 
applied  to  vessel  burst.s  given  data  following  testing. 

3.8  Evaluation 


Studies  into  fragmentation  began  in  earnest  in  the  1940' s,  centering 
almost  exclusively  on  the  fragmentation  hazard  presented  by  munitions 
blasts.  Much  classified  work  was  perl„..Tned  to  analyze  the  hazard 
presented  by  shrapnel  and  debris  to  personnel.  These  early  studies  are  of 
very  limited  usefulness  to  this  current  study  which  analyzes  the  effects 
of  fragments  produced  by  the  failure  of  a  pressurized  vessel. 

A  munitions  blast  tends  to  produce  many  more  sriller  fragments  than  a 
vessel  failure  due  to  overpressure.  The  increas'  in  the  number  o.*' 
fragments  produced,  however,  may  not  have  that  great  <au  impact  on  the 
maximum  velocity  achieved  by  any  fragment. 

There  has  also  been  a  fair  amount  of  research  into  the  failure  of 
vessels  containing  liquid  propellant.  This  is  of  primary  concern  in  the 
case  of  a  rocket  on  a  launch  pad.  Fragments  generated  by  the  part.iaL 
detonation  of  the  propellant  pose  a  hazard  to  ground  support  equipment  and 
may  also  lead  to  the  failure  of  nearby  liquid  propellant  vessels. 

Literature  on  a<_tual  failures  of  liquid  propellant  vessels  was 
researched  to  compare  their  fragmentation  with  that  of  p-ressur  i  zed  gas 
vessels  and  also  munitions  blasts.  Liquid  propellant  vessels  tend  to 
produce  fewer,  but  larger,  fragments  than  munitions  blasts  and  more,  but 
smaller,  fragments  than  gas  pressurized  vessel  bursts.  Only  the 
pressurized  gas  vessel  burst  is  considered  In  this  review. 

The  Taylor-Price  analysis  is  the  basis  for  most  models  which 
determine  the  velocity  of  fragments  generated  by  a  pressure  vessel  burst . 
Several  important  models  and  computer  codes  have  been  presented  which  have 
modified  the  Taylor-Price  analysis  to  improve  '..ne  assumptio  n.s,  or  have 
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uppMed  assumptions  to  analyze  a  special  case,  or  have  generalized  the 
assumptions  to  make  the  analysis  more  universal.  These  models  and  codes 
cover  a  wide  range  of  useful  applications.  The  biggest  drawback  is 

limited  data  to  verify  their  accuracie.s . 

There  has  been  very  little  data  published  on  fragmentation 
measurements.  Currently,  research  is  being  performed  by  Baum  to  measure 
fragment  velocities  and  derive  equations  for  worst  case  fragment 
ve.locit’es  for  a  large  variety  of  different  vessel  burst  cases.  It  is 
expected  that  he  will  publish  data  in  the  future. 

The  most  widely  used  fragmentation  data,  currently,  is  that  of 

Pittman.  However,  there  are  drawbacks  to  both  of  his  sets  of  data. 

In  the  1372  bursts  of  5  titanium  vessels,  only  2  vessel  bursts  were 
measured  by  both  the  breakwire  system  and  the  stroboscopic  cameras.  In 
one  of  these  2  bursts,  the  velocity  measurement  systems  differed  by  about 
25%  on  fragment  velocity.  Additionally,  Pittman  concluded  that  the  method 
for  estimating  the  number  of  fragments  produced  was  incorrect .  One 
important  conclusion  is  that  fragments  are  not  propelled  from  the  burst 
uniformly  in  all  directions. 

'.'he  results  of  the  197C  bursts  of  7  vessels  made  of  T-1  steel 

surprised  Pittman.  It  was  anticipated  that  the  ve.ssels  would  fail  in  a 
brittle  manner  and  generate  a  large  number  of  fragments.  The  fragment 
recr /ery  and  velocity  measurement  system  was  designed  accordingly. 

Instead,  each  vessel  failed  In  a  ductile  manner  and  burst  into  2 
fragments.  Four  of  t^o  14  fragments  were  not  recovered  and  were  a.ssumed 
to  have  traveled  many  hundreds  of  feet,  leaving  the  area.  Only  in  the 
Ccise  of  1  of  the  14  fragments  was  a  velocity  measured.  In  all  other 
cases,  velocitv  was  deduced  by  analyzing  the  pressure  -  temperature 
history  recorded  for  the  blast  wave  analysis. 

Another  matter  raised  by  tnis  experiment  is  the  question  of  how  a 
given  vessel  will  fragment  depending  upon  its  fail'..re  mechanism.  The 
failure  mechanism  could  very  well  play  a  key  role  in  determining  the 
number,  size,  mass  and  velo<;ity  of  fragments  produced.  To  date,  there  has 
not  beon  any  research  into  this  area. 
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Once  a  velocity  has  been  determined  tor  a  fragment,  calculation  of 
its  range  is  not  particularly  difficult.  The  ballistics  equations  have 
been  used  in  a  wide  variety  of  applications  with  success.  The  most 
critical  task  in  determining  the  range  is  determining  the  drag 
coefficient.  Drag  coefficients  have  been  measured  empirically  and 
tabulated  for  a  variety  of  objects.  The  choice  of  the  drag  coefficient 
presupposes  that  the  fragment  size  and  shape  are  known. 

Determining  the  distribution  of  fragments  produced  by  a  vessel  burst 
involves  cataloging  the  size,  shape,  mass,  speed  and  direction  of  travel 
of  fragments.  As  has  been  discussed,  there  is  no  model  or  theory  to 
predict  this;  studies  have  centered  on  statistically  analyzing 
experimental  results. 

The  only  usable  data  in  this  regard  is  the  1972  failures  of  5 
titanium  vessels  by  Pittman.  The  data  from  Pittman  is  not  sufficient  to 
draw  any  real  conclusions  beyond  the  fact  that  the  mass  di:  ‘  ribution 
appear.s  to  follow  a  log  normal  curve. 

More  experimental  failures  of  vessels  and  an  analysis  of  the 
mechanism  of  failure  will  be  needed  to  determine  fragment  distributions. 
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SBCTIOH  4 .  COKCLUSIOHS  AMD  RXCCMOCEMDATZOMS 


4 . 1  Summary 

As  described  in  this  report,  this  study  was  conducted  ic  (1) 
determine  and  summarize  existing  methodologies  for,  and  experimental 
confirmation  of,  energy  release  processes  (both  blast  and  f racjmentation) 
from  the  failure  of  pneumatic  pressure  vessels,  (2)  determine 
inconsistencies  in,  or  lack  of  information  on,  these  processes,  thus 
developing  requirements  for  future  study,  and  (3)  develop  a  plan  of  action 
to  resolve  inconsistencies,  or  generate  required  data,  to  establish  an 
effective  model  of  energy  release  from  failed  pneumatic  pressure  vessels. 
The  theoretical  approach  for  the  process  of  pressure  vessel  failure  has 
not  been  consistently  applied,  modelled,  or  experimentally  validated. 
Thi.s  results  in  a  significant  number  of  basic  questions  raised  as  a  result 
of  this  study  which  require  further  investigation.  It  also  results  in 
quest  10  13  raised  over  our  understanding  of  the  effects  of  failed  vessels 
in  aerospace  applications.  This  has  particular  importance  to  the  siting 
and  operation  of  vessels  on  permanent  manned  and  unmannned  orbiting 
platforms.  Therefore,  the  conclusions  and  recommendations  which  follow, 
although  originally  prompted  by  a  study  of  ground  support  equipment  (GSE) , 
have  direct  application  to  both  GSE  and  aerospace  vehicle  equipment  (AVE) . 
The  recommended  course  of  action,  if  implemented,  should  establish  a 
consistent  set  of  validated  models  for  energy  release  processes,  and 
subsequent  effect  on  surrounding  equipment,  for  the  failure  cf  pressurized 
vessels . 

4 . 2  Conclusions 

4.2.1  Blast  Wave  Methodologies 

As  presented  in  this  report,  any  of  several  blast  wave  methodol''.'gies 
(isentropic,  adiabatic  or  isothermal)  may  be  used  to  model  the  e.xpaosion 
of  gas  following  the  failure  of  a  vessel.  The  applicability  o  t  the 
methodologies  depends  on  the  falluie  environment  and  may  depend  on 
proximity  to  the  source  cf  the  failure;  for  example,  a  different 
methodology  may  apply  in  the  region  directly  adjacent  to  the  failure  (near 
field),  as  cc.-npared  to  a  regio..  more  distant  from  the  failure  (far  field)  . 
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Available  experimental  burst  data  indicates  that  current  standards  used  to 
predict  blast  wave  characteristics  may  be  overly  conservative  in  near 
field.  These  methods,  which  use  an  isothermal  expansion  model  coupled 
with  a  TNT  energy  equivalency,  may  also  be  underconservative  in  the  far 
field.  The  actual  expansion  of  a  vessel's  contents  following  rupture  is 
proposed  to  be  more  closely  approximated  by  an  adiabatic  process  rather 
than  an  isothermal  process.  However,  due  to  irreversibilities  in  the 
process,  the  use  of  the  adiabatic  expansion  equation  is  not  practical. 
Therefore,  an  isentropic  process  appears  to  be  the  most  applicable  model 
available  for  future  study.  Modifications  may  be  appropriate  as 

experimental  results  are  obtained  and  analyzed.  Evaluation  of  the  overall 
effect  of  using  ideal  or  real  gas  assumptions  must  al^o  be  included  in 
futur'-  study. 

Experimental  data  available  to  validate  existing  computer  models  and 
standard  approaches  is  generally  unavailable.  As  illustrated  in  Figure  4- 
1,  the  body  of  related  experimental  data  consists  of  11  burst  tests,  only 
two  of  which  were  made  under  an  8000  psi  burst  pressure.  In  addition,  no 
data  is  available  in  the  far  field.  Eval  ation  of  this  data  against 
available  models  does  not  yield  consistent  conclusiois . 

Two  computer  models  exist  Ic  simulate  the  r  y.  ;ansion  of  the  gas 
following  failure  of  the  ve.ssel  These  models,  titled  "Wondy"  and  "Wundy" 
(both  available  several  modified  versions) ,  although  used  to  predict 

blast  wave  character  atics,  have  not  been  validated  by  any  burst  tests. 

Based  on  th<!  limited  experimenta .  vessel  bu'  data,  the  follc  v.ing 
additional  items  remain  unaddressed,  tequii.ing  fuither  study  a.s  part  of 
the  recommendations. 

1.  Blast  wave  characteristics  -  Since  the  sh’tpe  of  the  o  'erpiessure 
wave  front  may  vary  .significantly  fiom  th.  predicted  by 
equiva lericv  wi'h  TNT,  an  evaluation  should  be  undertal^en  t  model 
and  mearure  shape  ^  hai  i:eris*'ic3  from  failed  pineu-  i.  v,  sure 

vessel.*^  . 
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PRESSURE  RANGES 

CYLINDRICAL 

SPHERICAL 

VESSEL 

VESSEL 

VERY  LOW  PRESSURE 

NO  DATA 

NO  DATA 

(<  350) 

LOW  PRESSURE 

640  PSI  (1.3) * 

NO  DATA 

(350  -  3000  PSD 

615  PSI  (1.6) 

MEDIUM  PRESSURE 

NO  DATA 

NO  DATA 

(3000  -  7000) 

HIGH  PRESSURE 

NO  DATA 

8145  PSI  (6.0) 

(7000  -  10,000) 

8015  PSI  (0.2) 

8015  PSI  (6.0) 

ULTRA  HIGH 

NO  DATA 

14, 965  PSI  (1.02) 

PRESSURE 

14,  765  PSI  (1 .02) 

(>10,000  PSD 

13. 415  PSI  (1.02) 

13, 815  PSI  (1.02) 

50. 415  PSI  (1.02) 

50,  415  PSI  (1  .02 ) 

*  (VOLUME  IN  CUBIC  FEET) 


Figure  4-1.  Pittman  Burst  Tests 
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2.  Burst  energy  component 3  -  Minimal  information  is  currently 

available  which  would  allow  evaluation  of  the  maKimum,  and 
most  probable,  percentage  of  internal  energy  available  to 
produce  the  blast  wave.  This  uncertainty  mny  yield  signi¬ 
ficant  variations  in  prediction  of  blast  wave  energies  and 
damage  potential.  Experimental  envelopes  should  be  developed 
for  distribution  of  energies  between  blast  wave,  fragmen¬ 
tation,  vessel  distortion,  and  crack  initiation  and 
propagation . 

3.  Jetting  direction  correction  factor  -  Since  nonuniform  blast  wave 

generation  is  most  probable  in  vessel  failure,  an  evaluation 
of  the  effect  of  jetting  on  the  items  addressed  above  should 
be  undertaken.  This  effect  may  result  in  a  significant 
increase  in  blast  wave  energies  &j.ong  one  radial  direction 
from  the  vessel.  Since  no  experimental  data  exi.gts  to 
address  this  effect,  data  should  be  taken  to  develop  a  model, 
practice  or  a  correction  factor. 

4.  TNT  conversion  factor  -  TNT  equivalency  is  currently  used  to 

relate  failed  pressure  vessel  blast  waves  to  blast  wave 
effects  on  equipment,  facilities  and  personnel.  Although  TNT 
blast  effects  are  well  defined  and  experimentally  confirmed, 
the  relationship  between  stored  energy  in  a  pneumatic  vessel 
and  an  equivalent  amount  of  TNT  has  not  been  experimentally 
confirmed.  In  fact,  the  literature  includes  a  number  of 
different  value,?  for  this  conversion.  Experimental  confirma¬ 
tion  of  the  applicability  of  this  relationship  should  be 
addressed  in  future  studies. 

In  summary,  there  is  a  significant  lack  of  complete  experimental 
data,  and  associated  validation  of  available  compute^’  mod"ls.  ir  appears 
that  current  prediction  methods  may  he  overconse  t  vat  i  vr>  in  come 
applications,  while  underconservative  in  other  applications.  ."jnce  the.se 
method.?  are  used  to  site  b.igh  pres.suie  v<-*.v.oel.s  and  sy.st  -it  gowr  rim'-.'it 

missile  and  space  complexes,  along  with  appl  i  ■  at  i  on.?  to  oil-itji.g 


74 


platfonoA,  it  is  prudent  to  address  these  issues  of  uncertainty  as  soon  as 
practical  and  resolve  discrepancies  between  current  standardized  methods 
and  actual  blast  effects. 

4.2.2  Fragmentation  Methodologies 

As  summarized  in  this  report,  the  study  of  fragmentation  mechanisms 
and  effects  have  been  studied  extensively  since  the  1940s,  with  initial 
emphasis  on  the  fragmentation  hazard  presented  by  the  detonation  of 
munitions.  The  failure  of  vessels  containing  liquid  propellants  has  also 
been  studied  extensively,  however,  not  until  the  early  1970s  have 
pressurized  gas  vessels  oeen  addressed.  The  Taylor-Price  analysis, 
presented  in  this  report,  is  the  basis  for  most  models  which  determine  the 
velocity  of  fragments  generated  by  a  pressura  vessel  burst.  This 
approach,  modified  to  address  a  variety  of  applications,  has  not  been 
satisfactorily  validated.  The  primary  conclusion  of  the  few  experimental 
rei,ults  available  is  the  fact  that  the  fragmentation  is  a  function  of  its 
failure  mechanism.  That  is,  the  failure  mechanism  may  play  a  significant 
role  in  determining  the  energy  imparted  (velocity) ,  number,  size/shape, 
and  mass  of  fragments  produced.  This  lack  of  information  creates  a 
significant  problem  in  using  these  models  as  a  predictive  tool.  The 
existing  models  are  useful,  however,  when  evaluating  failed  vesse-ls  when 
number  and  size  of  fragments  are  known.  Validation  of  models  should 
address  energy  distribution  between  blast  wave  and  fragmentation,  along 
with  ability  to  predict  resulting  effects  of  vessel  failure. 

4 . 3  Recommendat ions 

4.3.1  Short  Term  Recommendations 

As  described  above,  additional  data  is  required  from  actual  burst 
tests  to  validate  both  blast  wave  and  fragmentation  models.  This 
validation  would  include  the  m-idif ication  of  available  models  to  predict 
blast  wave  and  fragment  propagation,  thus  assisting  in  the  cibility  to  more 
accurately  locate  vessels  in  ground  support  and  orbiting  i/latform 
applications.  Short  term  goals  include: 
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•  obtain  additional  experimeiital  data  to  validate  and  assist  in 
modification  of  existing  computer  models  by: 

establishing  vessel  selection  criteria  for  burst  testing  a 
wide  variety  of  ground  and  flight  vessels 
inventorying  candidate  vessels  for  burst  testing 

•  develop  logistics  engineering  support  for  burst  testing  by: 

evaluating  test  site(8) 
identifying  test  support  requirements 
developing  cost  estimates  for  test  support 
developing  detailed  test  criteria 
developing  test  schedules 
developing  test  plans 

•  establish  baseline  computer  models  to  be  evaluated  for  blast  wave 
and  fragmentation  validation  by: 

obtaining  current  versions  of  blast  wave  and 
fragmentation  computer  programs 

running  computer  models  for  a  variety  of  initial  conditions 
to  establish  envelopes  for  experimental  parameters 

4.3.2  Long  Term  Recommendations 

Based  on  the  short  term  recommendations  described  above,  the 
following  long  term  reconanendations  have  been  developed: 

•  implement  a  series  of  burst  test  programs,  including  both  flight 
and  ground  support  vessels,  designed  to  validate  computer  models 

•  modify  computer  models  or  theoretical  bases  for  blast  wave  and 
fragmentation 

•  evaluate  relationship  between  TNT  equivalency  and  pressure  vessel 
burst  characteristics 

•  determine  validity  of  on-line  monitoring,  crack  propagation 
prediction,  fracture  mechanics  analysis,  blast  and  fragmentation 
effects,  as  applied  to  flight  a. id  ground  support  vessels. 
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4.3.3  Preliminary  Test  Plan  Matrix 


A  "strawman”  test  plan  matrix  was  developed  to  conceptualize  the 
considerations  necessary  to  implement  a  full  scale  vessel  burst  test 
program.  Figure  4-2,  tne  Preliminary  Test  Program  Matrix  for  Experimental 
Vessel  Burst  Assessments,  presents  an  initial  evaluation  of  burst 
configurations,  pressures  and  volumes,  along  with  a  variety  of  other  test 
parameters.  The  follov.'ing  paragraphs  describe  the  logic  behind  the 
selection  of  these  preliminary  test  programs.  Each  of  the  seven  test 
programs  described,  may  be  comprised  of  several  test  runs,  that  is, 
several  burst  tests  may  be  required  under  each  test  program.  In  addition, 
an  initial  trial  test  may  be  required  to  confirm  operation  of  the  test 
facility,  including  instrumentation. 

4. 3. 3.1  Initial  Burst  Test  Program  (Tests  1-4) 

Proposed  Test  Plans  1  through  4  are  designed  to  provide  a  baseline 
envelope  for  blast  wave  and  fragment  generation.  This  is  achieved  by 
controlling  (1)  the  failure  initiation  mechanism,  (2)  burst  pressure,  (3) 
number  of  fragments  produced,  and  (4)  symmetry  of  blast  wave  and 
fragmentation  pattern  produced.  This  control  is  a  result  of  the  use  of 
shaped  charges  to  initiate  failure,  removing  sufficient  wall  thickness  to 
allow  the  resulting  inadequate  remaining  wall  to  fail.  Each  of  these  four 
teat  plans  may  consist  of  several  individual  test  trials .  This  approach 
Can  also  be  utilized  to  (1)  validate  and  modify  existing  computer  models 
for  the  prediction  of  blast  wave  and  fragment  generation  and  propagation, 
(2)  develop  a  model  for  pred;ction  of  the  distribution  of  energy  between 
blast  wave,  fragmentation,  and  material  distortion,  and  (3)  modify  the 
test  facility  to  mor''  accurately  measure  t)>e  required  experimental 
parameters . 

As  illustrated  in  Figure  4-2,  Test  Plan  1  is  proposed  to  be  conducted 
using  flight  weight  vessels  with  a  volume  of  less  than  5  cubic  feet  and  a 
burst  pressure  of  approximately  1000  psi.  As  described  above,  since  tJiese 
tests  are  initiated  by  a  shaped  charge,  the  actual  design  pressure  for 
these  vessel.s  may  cover  a  wide  range  of  values.  The  primary  objective  ot 
Test  Plan  1  is  to  validate  and  modify  the  Wundy  blast  wave  computer  code 
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DESCRIPTION 


TEST  PLAN  1 


TEST  PLAN  2 


TEST  PLAN  3 
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1 

CoDput 
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TEST  PLAN  4 
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TEST  PLAN  7 


CYLINDER  (GSE) 
1000  psi 
25  cu  ft 


\  '  / 


1 

~n 

— 

1 

'  1  ' 


GN2 

Carbon  Steel 

Shaped  Charge 

Overpressure 
3-180O/90O 
Gas  Terapi^ress 
Frag  Vel,  Mass 

Blast , 

Wundy  Code 

TBD 

No 

No 
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TBD 

Yes 

Yes 

None  Defined 

Tests  2  &  4 

Also  Evaluate 
AET,t>l,  etc. 


CYLINDER  (GSE) 
5000  psi 
25  cu  ft 


\  I  / 


/■  'v 


GN2 

Carbon  Steel 

Machined  Flaw 

Overpressure 
3-1800/900 
Gas  Temp, Press 
Frag  Vel,  Mass 

Blast, 

Wundy  Code 

TBD 

Yes 

Yes 

None  Defined 

Tests  2,4  4  5 

Also  Evaluate 
AET,FM,  etc. 


CYLINDER  (GSE) 
10,000  psi 
25  cu  ft 


\  I  / 


/'\ 


GN2 

Carbon  Steel 

Machined  Flaw 

Overpressure 
3~1EOO/9CO 
Gas  Temp .Press 
Frag  Vel,  Mass 

Blast , 

Wundy  Code 

TBD 

Yes 

Yes 

None  Defined 

Tests  2,4,5  &  6 

Also  Evaluate 
AET.FII,  etc. 


Figure  4-2.  Preliminary  Test  Program 
Matrix  for  Experimental  Vessel  Burst 
Assessments 
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with  sscondary  objectives  of  (1)  validation  of  fragmentation  computer  codes 
for  the  simplest  two  fragment  case,  and  (2)  evaluation  and  modification,  if 
necessary,  of  the  test  array  orientation  and  configuration.  The 
circumferential  failure,  initiated  by  the  shaped  charge,  should  produce  a 
symmetrical  blast  wave  which  would  be  measured  in  a  plane  parallel  to  the 
failure  initiation  point.  This  configuration 

should  also  yield  the  maximum  blast  wave  energy,  with  the  minimum  energy 
inserted  to  the  two  fragments  and  with  minimum  material  distortion  prior  to 
failure.  This  should  be  used  to  begin  to  establish  the  upper  envelope  for 
the  blast  wave  energy. 

Test  Plan  2  introduces  failures  in  ground  support  pressure  vessels  utider 
the  same  conditions  as  described  in  Test  Plan  1.  The  blast  -j.’ /e  is 
reoriented  to  a  plane  perpendicular  to  the  horizontal,  with  the  two  fragments 
propagating  at  90®  to  the  blast  wave.  The  results  should  again  be  used  to 
validate  and  modify  the  respective  computer  codes  and  adjust  test  facility 
configuration . 

Test  Plan  3  is  proposed  to  use  flight  weight  vessels  at  a  burst  pressure 
of  5C00  psi  with  failure  plane  oriented  perpendicular  to  the  horizontal.  As 
described  in  Test  Plans  1  and  2,  two  fragments  would  be  produced  from  the 
failure.  This  program  should  complete  preliminary  evaluation  of  two 
fragment,  symmetric  blast  wave  failures  and  establish  baseline  envelopes  for 
the  remaining  plans . 

Test  Plan  4  initiates  multiple  fragmentation  using  shaped  charges  to 
cause  failure.  These  tests,  performed  on  ground  support  (or  flight)  vessels, 
should  be  used  to  estimate  the  partition  of  energy  between  blast  wave  and 
fragmentation  when  multiple  fraginents  are  involved.  Validation  of  computer 
models  should  continue  with  expansion  to  multiple  fragment  models.  This  test 
plan  is  designed  to  directly  compare  to  Test  Plan  2  establisliing  the  related 
envelope  which  should  be  the  basis  of  future  test  comparisons. 

Test  Plans  1  through  4  should  establish  the  validity  of  existing 
computer  models  and  provide  needed  data  for  modification  of  these  model.s. 
These  tests  also  establish  an  upper  and  lower  bound  envelope  tor  the 


partition  of  energy  between  blast  wave  and  fragmentation.  in  addition, 
re.3ults  of  the  tests  are  used  to  reconfigure  the  test  facility  and 
instrumentation  for  future  burst  tests  initiated  by  machined  flaws. 

4. 3. 3. 2  Immediate  Follow-on  Burst  Test  Program  (Tests  5-7) 

Proposed  Test  Plans  5  through  7  are  designed  to  .simulate  actual 
vessel  failures  using  machined  flaws.  These  plans,  with  burst  pressures 
at  1000,  5000  and  10,000  psi,  are  designed  to  address  several  significant 
issues  including: 

•  propagation  of  flaw  during  failure  and  ability  to  predict  fragmen¬ 
tation  pattern,  number  and  size 

•  validity  of  fracture  mechanics  in  predicting  flaw  growth  and 
approach  co  failure 

•  ability  of  acoustic  emissions  testing  to  predict  flaw  growth  and 
onset  of  failure 

•  validity  of  blast  wave  and  fragmentation  computer  codes  under 
actual  failure  conditions 

•  validity  of  energy  partition  envelope  for  blast  wave  and  fragmen¬ 
tation 

•  blast  wave  effects  on  code  predictions  including  jetting,  unsym- 
metric  blast  wave  distributions,  location  of  barriers,  etc. 

•  fragmentation  effects 

Details  for  Test  Plans  5  through  7  would  be  developed  as  part  of  the 
initial  burst  test  program,  and  as  such  only  major  goals  or  objectives  are 
described  here. 
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4. 3. 3. 3  Long  Term  Follow-on  Burst  Tests  (Test  8-  end  of  program) 

Specific  unresolved  issues  may  be  raised  as  part  of  the  initial  or 
immediate  follow-on  tests  which  should  be  addressed  in  longer  term  programs. 
These  would  be  developed  on  a  case-by-case  basis  and  may  be  Incorporated,  as 
appropriate,  into  existing  programs,  in  addition,  any  of  the  test  programs 
described  may  be  reordered,  modified,  deleted  and  replaced,  or  adjusted  due 
to  the  results  of  prior  tests. 

4.3.4  Summary  of  Recommendations 

A  test  program  should  be  initiated  to  address  several  major  is&  es 
raised  as  part  of  this  study.  These  issues  have  wide  appliC2d>ility  to  the 
location  of  vessels  on  manned  or  unmanned  space  platforms  and  as  ground 
support  equipment .  The  proposed  "strawman"  Test  Program  should  be  used  as 
a  basis  for  discussion  of  the  implementation  of  this  prograun.  Numerous 
related  research  efforts  may  be  incorporated  into  this  program  such  as  on¬ 
line  monitoring  in  space  and  ground  environments,  failure  effects  and  barrier 
design  options.  This  preliminary  evaluation  should  establish  a  point  of 
discussion  for  future  coordinated  full  scale  vessel  failure  prediction. 
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APPENDIX  A 


DEVELOPMENT  OF  TNT  EQUIVALENCE  FROM  DATA 


A-1 


Table  1.  Pittman  Data  from  {10  &  11] 


VESSEL 

NUMBER 

VESSEL 

SHAPE 

VESSEL 

VOLUME  (FT’) 

VESSEL  BURST 

PRESSURE  (FSI) 

VESSEL 

MATERIAL 

VESSEL 

MEDIUM 

1 

Cylindrica 

1  1.34 

640  T 

i  6-Al  4-V  Allc 

>y  GN2 

2 

Cylindrical  1.68 

615  Ti  6-Al  4-V  Alloy  GN2 

3 

Spherical 

0.238 

8015  Ti  6-Al  4-V  Alloy  GN2 

4 

Spherical 

6 

8015  Ti  6-Al  4-V  Alloy  GN2 

5 

Spherical 

6 

8145  T 

i  6-Al  4-V  Allc 

>y  GN2 

6 

Spherical 

1.02 

14965 

T-1  Steel 

GAr 

7 

Spherical 

1.02 

14765 

T-1  Steel 

GAr 

8 

Spherical 

1.02 

34415 

T-1  Steel 

GAr 

9 

Spherical 

1.02 

31815 

T-1  Steel 

GAr 

10 

Spherical 

1.02 

50415 

T-1  Steel 

GAr 

11 

Spherical 

1.02 

50415 

T-1  Steel 

GAr 

-2 


Table  Z  Expansion  Energy  •—  Ideal  Qas 


VESSEL 

NUMBER 

ISOTHERMAL 

ENERGY 

(FT-LB/FT’) 

ISENTROPIC 

ENERGY 

(FT-LB/FT') 

1 

0. 

35 

X 

10‘ 

0. 

15 

X 

10‘ 

2 

0. 

33 

X 

10‘ 

0. 

14 

X 

10‘ 

3 

7. 

25 

X 

10‘ 

2. 

41 

X 

10‘ 

4 

7. 

25 

X 

10‘ 

2. 

41 

X 

10‘ 

5 

7. 

39 

X 

10‘ 

2. 

45 

X 

10‘ 

6 

14 

.9 

X 

10‘ 

3. 

01 

X 

10‘ 

7 

14 

.6 

X 

10‘ 

2. 

97 

X 

10‘ 

8 

38 

.3 

X 

10‘ 

7. 

06 

X 

10‘ 

9 

35 

.1 

X 

10‘ 

6. 

5 

X 

10‘ 

10 

58 

.  9 

X 

10‘ 

10. 

42 

X 

10‘ 

11 

58 

.9 

X 

10‘ 

10. 

42 

X 

10‘ 

-3 


Table  3.  TNT  Equtvalencei  -  ktoal  Qaa 


VESSEL 

NUMBER 

TNT  EQUIVALENCE 

FROM  ISOTHERMAL 

(pounds  'INT) 

TNT  EQUIVALENCE 

FROM  ISENTROPIC 

(pounds  TNT) 

TNT  EQUIVALENCE 

FROM  PITTMAN 
l 

(pounds  TNT) 

1 

0.300 

0.132 

0.138 

(0.127) 

2 

0.358 

0.158 

0.169 

(0.156) 

3 

1.104 

0.366 

0.324 

(0.299) 

4 

28.17 

9.352 

10.1 

(9.33) 

5 

28.70 

9.510 

10.1 

(9.33) 

6 

9.83 

1.95 

0.82 

(0.757) 

7 

9.26 

1.92 

0.82 

(0.757) 

8 

25.28 

4.56 

0.95 

(0.878) 

9 

23.17 

4.21 

0.95 

(0.878) 

10 

38.98 

6.74 

1.0 

(0.924) 

11 

38.88 

6.74 

1.0 

(0.924) 

1  -  Based  on  1.545  x  10*  ft-lb,  per  pound  TNT  from  Kinney 


?.  -  Pittman  data  assumed  1018  Cal  per  gram  TNT.  The  values 
in  parentheses  assume  1.545  x  10*  ft-lb, . 


A-4 


Table  4.  Peak  Overpressure  —  Ideal  Gas 


VESSEL 

NUMBER 


DISTANCE 

FROM  BURST 

(ft) 

OVERPRESSURE 

FROM  ISOTHERMAL 

(psi) 

OVERPRESSURE 

FROM  TSENTROPIC 

(psi) 

3.3 

44.10 

23.90 

5.5 

14.37 

8.458 

12.0 

3.534 

2.386 

3.3 

50.21 

27.228 

5.5 

16.23 

9.44 

12.0 

3.87 

2.587 

3.3 

116.98 

50.96 

b.5 

36.67 

16.47 

12.0 

7.31 

3.90 

7.3 

213.5 

95.166 

12.0 

70.15 

31.49 

18.0 

28.62 

13.07 

7.3 

216.1 

96.66 

12.0 

71.22 

31.94 

18.0 

28.47 

13.226 

1.0 

2763.6 

1363.0 

10.0 

48.54 

14.9 

16.5 

16.5 

5.75 

60.0 

1.849 

1 

0.93 

Table  4.  Peak  Oyerpreaaufe  —  Meal  Qaa  (Continued) 


VESSEL 

NUMBER 

DISTANCE 

FROM  BURST 

<ft) 

OVERPRESSURE 

FROM  ISOTHERMAL 

(psi) 

OVERPRESSURE 

FROM  ISENTROPIC 

(psl) 

7 

1.0 

2700.42 

1225.6 

10.0 

46.43 

14.8 

16.5 

15.37 

5.7 

60.0 

1.804 

0.93 

9 

1.0 

3984.24 

1997.1 

10.0 

98.26 

27.27 

16.5 

31.58 

9.6 

60.0 

2.856 

1.3 

9 

1.0 

3852.6 

1931.0 

10.0 

91.96 

25.8 

16.5 

29.61 

9.2 

60.0 

2.73 

1.28 

10 

1.0 

4643.5 

2377.1 

10.0 

135.03 

36.5 

15.5 

43.42 

12.4 

60.0 

3.03 

1.57 

11 

1.0 

4  64  3 . 5 

2377.1 

10.0 

135.03 

36.5 

16.5 

43. 12 

12.4 

60.0 

,  ,  .i 

3.03 

1.57 

A-6 


Tabto  5.  Expansion  Energy  (H-lb/ft^  —  Real  Qaa 


VESSEL 

ISOTHERMAL 

ISENTROPIC 

NUMBER 

ENERGY 

ENERGY 

4.9  X  10* 
4.9  X  10* 
4.9  X  10* 
8.7  X  10* 
8.6  X  10* 


1.65  X  10 
1.65  X  10 
1.64  X  10 
1.76  X  10* 
1.76  X  10* 


TNT  EQUIVALENCE  (pounds  TNT) 

0.76  0.25 

19.26  6.41 

19.26  6.38 

5.75  1.16 

5.7  1,16 

PEAK  OVERPRESSURE  (psi) 

88.17  38.6 

27,63  12.67 

5.82  3,24 

717.1  395  9 

299.9  135.9 

52.7  23.3 

1^1.1  395.9 

299.9  135.9 

52.7  23.3 

2208.8  926.85 

32.4  8.82 

11.17  3.69 

1.47  0.66 

2208.8  926.85 

32.4  8.82 

11.17  3.69 

1.47  0.66 


APPENDIX  B 


TABLE  A- I 

from  Kingery  (9) 


B-1 


A 

AV 

X 

X 

AV 

2000 

6952. 

4 

7500 

1 

17  1  n 

2 

1800 

3 

1523 

2500 

5599 

4 

8000 

1 

1501 

2 

1900 

3 

1409 

3000 

4624 

4 

6500 

1 

1323 

2 

2000 

3 

1314 

3a00 

3097 

4 

9000 

1 

1182 

2 

2200 

3 

114K 

<►000 

3  34  1 

4 

9500 

1 

105') 

2 

2400 

3 

1U16 

4500 

2904 

4 

10.00 

9615 

1 

2600 

3 

9079-01 

5000 

2553 

4 

1100 

2 

8029 

1 

2800 

3 

8186-01 

5500 

2264 

4 

1200 

2 

6825 

1 

3000 

3 

74  30-01 

6000 

2022 

4 

1300 

2 

5920 

1 

3250 

3 

6640-01 

6500 

1818 

4 

1400 

2 

5186 

1 

3500 

3 

5980-01 

7000 

1645 

4 

1500 

2 

4665 

1 

3750 

3 

5430-01 

7500 

149  / 

4 

1600 

2 

4177 

1 

4000 

3 

4960-01 

8000 

I36fl 

4 

1700 

2 

3797 

1 

4500 

3 

4200-01 

8500 

1255 

4 

1800 

2 

3488 

1 

5000 

3 

3620-01 

9000 

1157 

4 

1900 

2 

3208 

1 

5500 

3 

3170-01 

9500 

1070 

4 

2000 

2 

2984 

1 

6000 

3 

2800-01 

UOOO  1 

99  33 

3 

2200 

2 

2596 

1 

6500 

3 

2500-01 

1100  1 

8602 

3 

2400 

2 

2299 

1 

7000 

3 

2260-01 

1200  1 

7544 

3 

2600 

2 

2061 

1 

7500 

3 

2050-01 

1300  1 

6670 

3 

2800 

2 

1867 

1 

6000 

3 

1B70  -01 

1400  1 

5923 

3 

3000 

2 

1706 

1 

9000 

3 

1580  -01 

1500  1 

5334 

3 

3250 

2 

1537 

1 

1000 

4 

1370-01 

1600  1 

4782 

3 

3500 

2 

1397 

1 

1700  1 

4322 

3 

3750 

2 

1279 

1 

1800  1 

3919 

3 

4000 

2 

1178 

1 

1900  1 

3540 

3 

4500 

2 

1015 

1 

2000  1 

3207 

3 

5000 

2 

8876 

2200  1 

2630 

3 

5500 

2 

7857 

2400  1 

2180 

3 

6000 

2 

7023 

2600  1 

1834 

3 

6500 

2 

6328 

2800  1 

1558 

3 

7000 

2 

5742 

3000  1 

1337 

3 

7500 

2 

5222 

3250  1 

1117 

3 

8000 

2 

4769 

3500  1 

9438 

2 

9000 

2 

4041 

3750  1 

6064 

2 

100.0 

3 

3484 

4000  1 

6958 

2 

1100 

3 

3047 

4500  1 

5316 

2 

1200 

3 

2692 

5000  1 

4164 

2 

1300 

3 

2405 

5500  1 

33  76 

2 

1400 

3 

2162 

6000  1 

2782 

2 

1500 

3 

1970 

6500  1 

2334 

2 

1600 

3 

1793 

7000  1 

19,89 

2 

1700 

3 

1047 

B-2 


PiMman  Vessel  1 
Ideal  QN2 


C-2 


o 


Pittman  Vessel  3 
Real  QN2 


C-5 


(psi) 


Pittman  Vassal  4 
Ideal  QN2 


C-6 


Pittman  Vaasel  4 
Real  GN2 


DISTANCE  (It) 


10 


OVERPRESSURE  (psi) 


Pittman  Vessel  5 
Ideal  QN2 


4 


Pittman  Vessel  5 
Real  QN2 


c.» 


100 


(psi) 


I 


Pittmitf)  Vms«I  7 
Ideal  QAr 


OVERPRESSURE  <psi) 


OVERPRESSURE  <psi) 


Pittman  VasMi  8 
IdMl  QAr 


C-12 


(psi) 


Pittman  Vessel  9 
ideal  QAr 


OVERPRESSURE  (psi) 


10,000 


Pittman  Vassal  10 
Ideal  QAr 


C-14 


OVERPRESSURE  (psi) 


10,000 


Pittman  Vessel  11 
Ideal  GAr 


f 


C-15 


